


JOURNAL OF THE 
AERONAUTICAL SCIENCES 





VOLUME 10 


FEBRUARY, 1943 


NUMBER 2 





Three-Dimensional Wing Flutter Analysis 


ALEXANDER H. FLAX* 
Curtis-Wright Corporation 


ABSTRACT 


A set of equations is derived for the flutter of a wing-aileron 
system in three dimensions using methods familiar to the average 
The aerodynamic force coefficients are obtained by 
The results are 


engineer. 
strip integration of two-dimensional coefficients. 
presented in a form suitable for direct calculation by assuming 
that the fundamental vibration modes of the wing in bending and 
torsion may be used to represent the corresponding displacements 
of the wing in the flutter motion. The validity of the assumption 
made is discussed. 

Empirical weighting methods commonly used in two-dimen- 
sional analyses to account for three-dimensional effects are con- 
sidered in the light of the results obtained in the paper. Approxi- 
mate equations are given for the correction of coupled frequencies 
obtained in still air vibration tests to the “uncoupled” frequencies 
desired for flutter analysis. The corrections for the effects of 
the apparent inertia of the air are also included. 


INTRODUCTION 


T ‘HE PRESENT TRENDS toward higher speeds and 
lighter structures have made flutter considerations 
increasingly critical in the design of structural com- 


three-dimensional 
two-dimensional 


ponents. The reduction of the 

wing structure to an analogous 
structure, which has in the past been the most com- 
monly used method of analysis, has proved to be un- 
satisfactory for many practical problems. Local 
variations in weight, moment of inertia, and c.g. loca- 
tion have made the choice of suitable two-dimensional 
parameters extremely questionable. Large local vari- 
ations of this type frequently arise in wings that are 
used to carry fuel or armament. Arbitrary methods 
of weighting parameters have also proved objection- 
able, because these methods have often been based on 
unsound assumptions (as will be shown later). 

A completely satisfactory solution of the three- 
dimensional flutter problem can be obtained only by 
the development of a dynamically and aerodynamically 
sound theory for flutter in three dimensions. Un- 
fortunately, no complete solution of the aerodynamic 
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problem for a vibrating wing of finite aspect ratio is 
yet available. A partial solution (no aileron forces 
are considered) is given in reference 9, but it seems to 
be too cumbersome for adaptation to practical analy- 
sis. Thus, at present, all consideration of aerody- 
namic forces on the wing must be based on strip 
theory—that is, the forces on each spanwise element of 
wing are computed by assuming that the flow there is 
independent of the fluid motion at all other points on 
the wing. This assumption involves neglecting the 
induced velocities caused by the free and trailing 
vortices all along the span. 

Lyon® states that the effects of these induced veloc- 
ities is not as great in a periodic motion of the wing as 
in a steady motion. This is attributed to the fact that 
in a periodic motion the successive free and trailing 
vortices are of opposite sign, and their induced veloc- 
ities tend to cancel one another out except in the 
immediate neighborhood of any given vortex. Thus, 
in a vibrating wing, the effect of finite span will be due 
mainly to vortices immediately behind the wing. 

The most serious errors occur near the tip, as is 
pointed out by Cicala’. This is to be expected, be- 
cause not only are the trailing vortices largest near 
the tip but also, in a wing vibrating near the funda- 
mental mode, the amplitudes, which determine the 
aerodynamic forces, are the largest near the tip. The 
induced velocities tend to reduce the lift on the wing 
and are equivalent to an aerodynamic damping. The 
effect is to increase the flutter speed so that neglect of 
span effects is conservative. 

For wings of large aspect ratio, neglect of aero- 
dynamic span effects is not likely to lead to large errors 
in calculated flutter speeds as was pointed out by 
Theodorsen and Garrick’. For tail surfaces, however, 
where the aspect ratio may be three or lower, the 
error may be quite large. This is borne out by the 
fact that the flutter modes involving tail surface bend- 
ing torsion have led to flutter speeds, when calculated 
by conventional methods that had proved to be success- 
ful on wings, which were far lower than practical ex- 
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perience would admit. For tail surfaces, it may be 
possible to apply empirical corrections to the aero- 
dynamic force coefficients to account for aspect ratio 
effects. Such phases of the problem as have already 
been solved (presented in references 8 and 9) should be 
of considerable value in determining the empirical 
corrections. The equations set up in this paper allow 
for inclusion of any such corrections to the local aero- 
dynamic force coefficients. 

A rigorous solution of the problem of three-dimen- 
sional flutter insofar as the dynamic and elastic phases 
are concerned offers no fundamental theoretic diffi- 
culties, since the methods of classic mechanics are 
sufficient to solve the problem. The difficulties in- 
volved are purely those of practical computation. 
The wing in three dimensions constitutes a system 
with an infinite number of degrees of freedom. [If all 
wing displacements are assumed to consist of bending 
displacements, /, of the elastic axis, and torsional dis- 
placements, a, about the elastic axis, then the defor- 
mation of the wing may be described as a function of 
time and spanwise location by the two infinite series 
(if a solution that is a sinusoidal function of time is 


assumed): 
h(x, t) = aye f,(x) + are + fox) +2... 
a(x, t) = dye**™ gic) + doe *™ go(x) 22... 


Where the a's and b’s are constants, the f’s and g’s are 
arbitrary functions for the shape of the displacement 
curves that satisfy the boundary conditions in bending 
and torsion, respectively, and the 6’s and ¢’s represent 
the time phase relationships between the various terms; 
w represents the natural frequency of the periodic 
motion. Each term in the series represents a degree of 
freedom, and the accuracy of any work depends on the 
number of terms employed. The wing motion could 
probably be sufficiently well described by the use of 
two or three terms in each series, but when at least one 
degree of freedom is added to describe aileron motion, 
a total of from five to seven degrees of freedom results. 
Inasmuch as solution of the flutter problem involves a 
trial-and-error solution of a determinant of order equal 
to the number of degrees of freedom, such a procedure 
is seen to be unsuitable for routine analysis. 

An alternative method proposed by Lombard,' as 
well as by various British investigators, is to use only 
two degrees of freedom to describe the displacement 
of the wing, one for bending and one for torsion. The 
displacement curves are assumed to be the same as 
those obtained from vibration of the wing in the 
fundamental modes for bending and torsion. These 
curves may be obtained either from ground vibration 
tests or by calculation.'' This procedure implies that 
flutter occurs in the fundamental mode. This, of 
course, is not necessarily true, but for wings having a 
damping coefficient comparable to that of most modern 
metal wings, the error will not be serious as pointed out 
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by Theodorsen and Garrick.’ In this paper, the 
fundamental modes of bending and torsion will be 
assumed to be the displacement curves in the flutter 
motion, in order that the number of degrees of freedom 
may be limited to three. 

Several papers on three-dimensional flutter analysis 
have already been published; those of Lombard,! 
Loring,” and Bleakney* are of particular interest. 
Lombard has derived a set of three-dimensional equa- 
tions using assumed deflection curves. He has not, 
however, employed the commonly used notation of 
references 6 and 7 and, further, his analysis does not 
include aileron effects. Loring and Bleakney have 
both derived sets of general equations applicable to any 
number of degrees of freedom which undoubtedly con- 
stitute important contributions to the basic theory of 
flutter in three dimensions. Loring, however, has em- 
ployed instruments of analysis, such as Lagrange’s 
equations and matrices, which can hardly be regarded 
as part of the equipment of the average aeronautical 
engineer. Bleakney too has employed a method of 
derivation not familiar to the average engineer. The 
results of both these papers have not been developed 
to the point where immediate use for routine calcula- 
tions is possible, mainly because both authors wished to 
allow for treatment of a large variety of problems in- 
volving many possible degrees of freedom. 

It is the purpose of this paper to develop, by using 
principles of mechanics and mathematics within the 
scope of the average engineer, a method for three-dimen- 
sional wing flutter analysis in three degrees of freedom. 
The nomenclature of references 6 and 7 is employed 
wherever possible. 


THE EQUATIONS OF MOTION 


Consider a spanwise element of wing and the forces 
acting upon it as in Fig. 1. The case of torsional 
motion is chosen as typical. If a is the angle of twist, 


x 







CH +8 (C# dy 


FLASTIC_AXIS 














Fic’ 1. 


then the torque on each side of the element is COa/Ox, 
where C is the torsional stiffness. In addition, there 
are body forces in the element equal to the inertia con- 
stants times the corresponding accelerations, as well as 
surface aerodynamic forces proportional to the dis- 
placements, velocities, and accelerations of the ele- 
ment. The body inertia forces (per unit length) are 
given on page 9 of reference 6 as: 
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02a oh where a’, 8’, and h’ are functions of spanwise distance 

I, + Up + We — 2) sis : e+ tore P 
or? ot? from the wing root, x, alone and represent the temporal 
maxima of the displacements. ¢; and ¢2 are the phase 


Assuming simple harmonic motion for the displace- angles by which the aileron and bending displacements, 


ments, respectively, lead the torsional motion. Then 
a= a'e™ 0?a/0l? = — wa 
ad (wt + d1) 2) 2 = 2 
8 as B’e i(wt 1 078/Odt = —wWw B 
h = hietor+#) 0°n/dt? = — wh 


Summing all moments about the elastic axis of the element, and noting that since time has been eliminated as a 
variable the partial derivatives may be replaced by total derivatives leads to 


— ¢(c%) — wI,a — w*[I, + b(c — a)S3]8 — w?S,h — M, = 0 (1) 
dx dx 


where /, is the moment pic the elastic axis of the aerodynamic forces. This will be treated later. Maultiply- 
ing through Eq. (1) by a/2 and integrating with respect to x from 0 to /, where 0 is at the root of the wing and 
l at the tip, gives 


1 1 2 1 1 
- ff: s (c 42) ax - wf L—&- of Usp + b(c -— a)) 59 < dx _ of S. Fe -f M,< dx =0 
0 2dx\ dx 0 2 0 2 0 2 
(2) 


The first term of Eq. (2) may be further reduced by partial integration using the formula 


[ru B ac = lwo — fo ax 
x dx a dx 


Let a/2 = uand tA (c%2) = dy 


dx dx, dx 
1 
Then f = £ (c dc) dx = 
0 2 dx dx 


t 2 
ac fel - [&(%) dx 
2 dx \o 0 2\dx 


Since a = 0 at x = 0 and da/dx = 0 at x = I, the first term on the right side of the equation above drops out, 


giving: 
1 1 . 
f pd (c%) en ~ £ : (42) dx (3) 
0 2dx\ dx 0 2 \dx 
Substituting Eq. (3) into Eq. (2) 


1 9 1 2 7 i 
f £(#y dx — oft — de — wf Tg + d(c — a) Sp] > Bas — coat f g. he 3. — v4 M,< dx = Q 
dx 0 2 0 2 0 2 


(4) 


The expressions h = h,fi(x) and a = a;,fo(x) may now be substituted in Eq. (4), where f; and fe are the funda- 
I / ; 

mental vibration modes in bending and torsion, respectively, and the subscript ¢ indicates tip value of the dis- 

placements. Eq. (4) becomes after division by fo'(I,a?/2 )dx, 


l ia 1 fod x F 
Wot — wa, — wip Lo Ua + be = aS fide _ oy Se'Safifde _ JaSh'Maade _ 4 (4a) 


ST fo°dx " S, Ta fa2dx ST, fedex 











since the expression 
So'Cl(fo’)?/2) dx/ Sl '(Iafe?/2)dx = w,? 


where .w, is the fundamental torsional frequency involving the f, mode (see page 185, reference 5). 
By similar treatment of the equations of equilibrium for aileron moments and for bending forces as given on 
page 9 of reference 6, two additional relationships are obtained. 
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The equation for aileron motion becomes 


2 “yal \2 yal 2 J “yal oes ia *yal 
r - Way Way w-h;y : 
(wg? — w?)B — — / I3fodx — “/f b(c — a)Sgfodx — | Safidx — / M,8dx = 0 
™ I yl. yu af I; yl yil a Igy yiul ois Igy B yil r 


(4b) 


where y/ is the length of the aileron and y,/ and y2l are the distances from the wing root to the inboard and out- 
board ends of the aileron, respectively. 
The bending equation leads to: 
: ’ on cle Salas o. de Sabsfrudx 1 Sf-'Phdx 
(w,? — w)h, — w?B ———=—_. — wa, —. ie - = = 0 (4c) 
J, 'Mfirdx JS, 'Mfrdx h, J,'Mfirdx 


AERODYNAMIC FORCES 


Theodorsen on page 9 of reference 6 presents the equation for the forces and moments on a unit strip of airfoil 
as linear functions of the displacements. 
These equations may be expressed in the following form (first suggested by Loring”): 


oe ee oe Oe ee 
M, = —pw'b?(A,,h + Aga ba + Agg 0B) (5) 
My = —xpw'b(Ay, h + A,, ba + Arg bf) J 


where A,, is equal to R, + t,; on page 12 of reference 6 with all mass and elastic terms omitted (terms con- 


taining K orQ). Then, 


I rf 
M,,adx Fd 1 ‘ 1 4 "yt /p\4 
3q)* t * » 120 ) 6 | . h e ( 
— . i. : k [ (°) AuSifadx + = . I (: ) A jafo*dx + Bbc / ( .) Aus fats 8) 
JS, T fo*dx JS, Te fe*dx be 0 bo l 0 bo l J vil bo 


0 


“yal ow t “ 3 . c Do sai ; . bo nied ; 
wa 4 | Mepis = fF / (°) hte 4 Se [ (*.) Ay. fad: + [ (°) Angi (7) 
Igyl B. yil Il l yil b, l J yil bo l J vil bo 


Finally 


[ 
Phdx . as 
h, 0 o"w* h, ‘ . ocak abo ’ 3 : Bho ve /b \3 2 : 
ee * J (°.) Aa fidx + 22 7 (*.) Awhifdx + © [ (°) Ag fits | (8) 
JS, 'Mfirdx JS, 'Mfirdx L Jo \bo l 0 \bo L Sa Nbo 


9 


It will be noted that in Eqs. (6), (7), and (8), an arbitrary reference semichord d) has been introduced. An exami- 
nation of the equations will disclose that the particular chord chosen has no effect whatever on the overall values 
of the functions on the right side of Eqs. (6), (7), and (8). The A’s for any given wing are functions of V/bw where 
b is the local semichord at any station. In solving for the flutter speed values of the aerodynamic force, integrals 
will be required for several values of V/bww. This necessitates recomputing the integrals several times graph- 
ically because the A’s are complicated functions of V/bw, which of course varies from station to station. The 
difficulties involved in such calculation may be circumvented by using simple approximate expressions for the 
A’s. These expressions may be obtained by fitting the curves of A versus V/bw for the given wing obtained from 
reference 6. The A value may be broken down to real and imaginary parts, so that A = A, + 7A,. It has 
been found that in the range from V/bw = 0 to V/bw = 7, the value of either the real or imaginary part of A 
may be satisfactorily approximated by (at most) a parabolic function. Thus, 


Ap = Ag, + Mp (V/bw) + ng(V/bw)? 


Substituting into an integral of the form: I, = f, '(b/bo)" Ag f(x)dx gives 


Gy Las +m ) + eZ) Joe 
I, = — Ar + mz | — r x)dx 
. P (> : ~ bw ™ bw Sex 
which may be written as 
l n , i n—l rT\2 I n 2 
I, = Ag f (5) f(x)dx + mz (“) f (°) f(x)dx + np (“) Z (;-) f(x)dx (9) 
0 bo bow 0 bo bow 0 bo 
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The integrals in Eq. (9) do not vary with V/bow and need only be computed once. Thus, the computation of 
the aerodynamic force terms for various values of V/bow, once the integrals are initially computed, is no more 
difficult than in the two-dimensional case. 

In some cases, especially where the variation of } along the span is not great, it may be possible to regard the 
A’s as constant over the entire span and to compute them on the basis of the value of V/bw at the three-quarter 
A simplification of this type seems especially justified because the A values are themselves 


semispan section. 
Lombard! makes an assumption of 


necessarily approximate, neglecting as they do the effects of finite span. 
this type in his equations. Further work on this matter is necessary with a view toward determining a simple 
yet close approximation, possibly even including span effects in a rough manner. 


SOLUTION OF THE STABILITY DETERMINANT 


Substituting Eqs. (6), (7), and (8) into Eqs. (4a), (4b), and (4c), the equations of motion are reduced to the 
form 
Aa, + BB + Ch, = 0 (torsion) 
Da, + EB + Fh, = 0 (aileron) (10) 
Ga, + H8 + Ih, = 0 (bending) 


II 


The coefficients of a,, 8, h, are complex functions of the physical parameters of the airplane, w,/w and V/bw. 
Damping coefficients g,, gz, and g, may be introduced in a manner similar to that used in reference 7. The fol- 


lowing symbols are employed in the coefficients: 
M, = JS, Mfdx I,’ = Igy ) = So Ta fe2dx 
Eqs. (4a), (4b) and (4c) are divided through by w?, and the coefficients become, 


2 3 l \4 
A= [2 + ign) — 1] + OH a eae | 
aw i Ll Jo bo 


Tr 


yl 
f [Ig + b(c — a) Sg] fodx ‘ —_ r 
B= yal + col f (>) Ausf 
Y 0 


i I . tise © 


Tr Tr 


l 
7 f Safi fedx mi aa f (; Seance 
0 


1 . 2s Y 


bos 


JO ts + We = a) Spi fete ; ibis 
a ee et / @zeza 


ar 2 L-Lit«u WV 


2 3 yal 4 
ht « E (1 + igs) — | Pe |* fi (*) Angi (10) 
w? I,' L Sw \bo 
| % , 
Sp fidx 3 yal 3 
oa me yil ae oe + mpbo l k J (? ) A, fidx | 
I,’ I,’ l yil bo 
l 
ite Js! ee mobi k if “€ J Aa fife) 
M, M, l Jo bo 
js 
Ss fidx » Fol 2 
a yil ai tpbo"l > f ie (¢ y 4 ax 
sli a. oe F »y 


2 blll ('/b\2 
I=|% t+ ies) — 1] + OE LEY aa peae 
ex “ M, Lido \d, pom 
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Eqs. (10) will all be satisfied only when the deter- 
minant of the coefficients is equal to zero or 


ABC 
DEF| = 0 (11) 
GHI | 
This leads to the relationship 
AEI + BFG + CDH — CEG — AFH — BDI = 
(12) 


The unknowns in Eq. (12) are V/b,w and w/w,; the 
real and imaginary parts of Eq. (12) must each in- 
dependently be equal to zero so that two independent 
equations are provided for determination of the two 
unknowns. The procedure followed is to assume several 
values of V/bo» and to solve for w/w, for each value 
from both the real and imaginary equations. Two 
curves of V/bow versus w/w, are thus obtained, one 
from the real part of Eq. (12) and the other from the 
imaginary part of Eq. (12). The values of (V/bw),; 
and (w/w,); at the intersections of the curves represent 
critical flutter speeds. The flutter speeds may be 
found as 

(13) 


V = (V/bw) (w/w) dow 


RELATION TO Two-DIMENSIONAL ANALYSIS 


It will be of some interest to compare the results of 
the three-dimensional analysis with some of the com- 
monly used arbitrary methods for putting three-di- 
mensional effects into the two-dimensional methods. 
Consider the commonly used weighting formulas as 
given in reference 10. 


So'Mhdy 


M = = (14) 
S'hdy 
where JM is the mass per unit length, and 
= Si pady/ fr! ady (14a) 


and p is the radius of gyration about the c.g. expressed 
in per cent of chord. 
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Examination of Eq. (4c) shows that ./ should be 
weighted as h? rather than as h. Weighting in propor- 
tion to # tends to overemphasize inboard parameters 
sometimes leading to appreciable errors. The same 
argument may be seen to hold for Eq. (14a) from 
Eq. (4a). Even if the weighting is carried out as the 
square of the deflection curves, discrepancies still re- 
main. For instance, the value of x, actually must be 
weighted differently for use in the different equations. 
This can be seen in Eqs. (4a), (4b), and (4c). Terms 
in Sg, appear as follows: 


InH/ S,X"Sgfidx/M, 
InB/ f¥"(c — a)bSgfedx/I, 
InF/ f,¥"b(c — a)Sgfidx/T,’ 


The two-dimensional theory allows for only one value 
of xg unless the analyst rederives the xg terms and re- 
vises the equations. 

If the weighted chord is used for the determination 
of aerodynamic properties, the problem of which power 
of the chord to use in weighting arises. Examination of 
aerodynamic terms in the equations of this paper shows 
variation as the square, cube, and fourth power of the 
chord in different terms. 

The equations of this paper, being three-dimensional, 
include the effects of partial-span ailerons. Theodorsen 
and Garrick’ modify the two-dimensional analysis for 
this effect by means of a coupling factor £, applied only 
to the aileron equation. The terms D and F indicate 
that for three degrees of freedom the terms of aileron 
equation must be modified by two coupling factors, one 
based on / and the other on a. In addition, the terms 
giving the aileron effect in the bending and torsion 
equations must be modified by similar coupling factors 
as in terms B and H. 

It is hoped that this brief discussion will aid in the 
intelligent choice of two-dimensional parameters for 
use in various preliminary analyses where a three- 
dimensional analysis is not justified because of sketchy 
data. 


GROUND VIBRATION TESTS 


The methods developed in this paper will be of considerable value in interpretation of ground and vibration 


tests. 


the fictitious “uncoupled” frequencies desired for use in the flutter equations. 
effects of apparent air mass which are not desired in the “uncoupled”’ frequencies. 


The tests always give frequencies that represent couplings between bending and torsional motion instead of 


In addition, the tests include the 
The solution for the ‘‘un- 


coupled” frequencies in terms of the coupled may be obtained from Eqs. (4a) and (4c) set down for V/bw = 0. 


Then Eq. (4c) gives 


a eae 


and Eq. (4a) gives 


Se Safefidx 
0 'M(fi)°dx 





We ? h, ry fifedx 
— lla, -—- ——=——— + 
( ) | ST fa)*dx 


w 





I 2 l 3 
0 0 
apbo'l Ae f( 4) ( ) ade ("(2) ‘ rs 
_ aes pee a - sj] — 9 2d4 aS 0 if a 
oi ba ST) Te 


1 pDo"l | 
CM ( fy)2dx 


(15) 
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The integrals will all be needed in any event for use in the flutter analysis so no additional evaluation of integrals 
is involved. Then the equations may be written: 


| (#)’ +Q | h, + Pa; = 
Mh, + | (**) 4+ N | a, 


The determinant of the coefficients of these equations must equal zero in order that they may be satisfied for all 


values of h, and a, This leads to the equation 


ae ae (=) | oe 7: a ee ve) ( ep ae 17 
(<) 1 (2: G+Ni = ¥[(“) +N] -4(* (ON PM) ¢ /2(QN PM) (17) 


Thus, two coupled frequencies, w; and w,, may be computed from the calculated ‘‘uncoupled’’ frequencies w, and 


To solve for the “‘uncoupled”’ frequencies 


w, and may be checked against coupled frequencies obtained in the tests. 
These are: 


from the coupled, #, and w, are substituted in Eq. (7) for w, giving two equations in , and w, . 


(2) (*) 4 (**) O+ (*) N+QN—PM=0 

Wy WW Wy ®) (18) 
(2) (*) $- (**) Q+ (2) N+QN — PM = 

We We We We 


Values of w, are assumed and w, is calculated from each equation. Two curves are plotted and the point of 


intersection represents the solution. 
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Rate of Climb Calculations 
G. OTTEN* 


Bandoeng Unwersity 


INTRODUCTION 


REDERICK C. PHILLIPS! has pointed out that in vari- 

ous cases the assumption of zero acceleration along 
a flight path is not correct and has given an equation 
by which the rate of climb calculations made in the 
usual way could be corrected for the error of neglecting 
the dynamic effect of the acceleration. 

It is also generally known that, especially for air- 
planes of great climbing capacity, the assumption that 
in the position for climb the lift should be equal to the 
gross weight causes other inaccuracies in the calcula- 
tions of the rate of climb.’ 

The correlations that have to be made if no account 
is taken of the acceleration or the smaller amount of 
the lift when the airplane climbs will generally be in an 
opposite direction and partially eliminate each other. 
While in such calculations the results are not unsatis- 
factory, decreasing tolerances, increasing engine powers, 
higher wing and propeller loadings by decreasing drag 
coefficients make it desirable to look for a better method 
of calculating the rate of climb.? 

With the method described below and making use of 
normal data, the climbing speed of every airplane can 
be determined with great accuracy without materially 
increasing the amount of work involved. 


SYMBOLS 


The following symbols will be used: 


F = force in general 

W gross weight (Ibs.) 

1 = thrust (Ibs. 

D drag (Ibs. 

p = air mass density 

f = airplane drag area 

q total energy pressure = pV?/2 

g = acceleration of gravity 

0 inclination of flight path to horizon (radians) 

C = flight rate of climb (ft. per sec 

V air speed 

t = time (sec.) 

A = horsepower equivalent (550 ft.lbs. per sec.) 

(hp = engine power 

n = propulsive efficiency 

C, = lift coefficient correctel for the whole airplane, 
when climbing = L/qf 

Cp = drag coefficient corrected for the whole airplane, when 
climbing = D/qf 

C,, = lift coefficient for horizontal flight = W/gqf 

Cp; = drag coefficient for horizontal flight 

C, = climbing speed coefficient = CW/qVf 
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DERIVATION OF THE FORMULAS 


The power delivered to the airplane by the propeller 
Ayn (hp.) can be divided into the power DV necessary 
to overcome the air resistance, the power WC necessary 
for climbing, and the power (VW/g)(dV/dt) for the 
the airplane, so that: 


= DV + WC + V (W/g)(dV/dt) 


acceleration of 


An(hp.) (1) 


and because C V sin 0 


An(hp.)/V = D+ W sin 6 + (W/g)(dV/dt) (2) 
As, according to Fig. 1, the forces acting perpendicu- 


larly to the direction of the flight path must be in 


equilibrium: 
L = W cos 0 (3) 
hence, 
W sin 6 = V (W? —'L?) (4) 
L Ary (bp) 
¥ V 











Fic. 1. 


Using Eq. (4) and the well-known aerodynamic equa- 
tions: 


D = Coqf 


and 
is = Cigf 
Eq. (2) may be written in the following form: 
An(hp. er W\ (dV ; : 
fee) — Coqf - ( )(' ) = V(W? — Cr7q*f*) 
J g dt / 


(7) 











RATE OF CLIMB 


By squaring both sides of Eq. (7) and rearranging: 


An(hp.) W\(/dV 
Ct + Co? — 2Cp | Asp. 5 ( )(' )| * 
Vof gqf/ \ dt 
| 428. ne (=\=)I  _ 
Vof gqf/\ dt qf 
Supposing that the airplane is flying at a certain 
selected altitude with the throttle entirely open and 
that the engine is giving its maximum power, then Eq. 
(8) only has terms that are dependent upon a variable 
V. 
Eq. (8) therefore represents a series of circles in the 
C,-Cp plane, the centers of which are on the Cp axis 


Further considering the circles, shift the C, axis in 
such a way that this axis goes through the points 


= 2p) _ (3 (4 r) and C, = 0 (9) 
Vof eqf/\ dt 


then, by arranging Eq. (8) with the help of ’ 


An(hp.) (W (e 
Vaf (=) =) 
the equation for the circles becomes 
Cr? + Co? — (W?/g*f?) = 0 
Further, assuming that 


W = Cin qf 


=0 (8) 


Cp 


Cp = Cor -+- 


(10) 


(11) 


in which C,, is the lift coefficient when flying hori- 


zontally at the speed |’, Eq. (10) will read as follows: 
G+ 

As C_, is the lift coefficient of the weight W of the 
airplane and C;, is the coefficient actually used and of 


Cor? = Cr’ (12) 


the lift L, Eq. (12) gives some idea of the errors that 
are made by assuming L = JW. 

A circle of the form of Eq. (12) described in the form 
of Eq. (8) will intersect the Cp-C, curve of the air- 
plane, the polar curve, in zero, one or more points. 

By not considering those intersections that yield a 
negative lift coefficient and those situated in the area 
in which the wing is stalling, then, with usable airfoils, 
only one or zero intersections will be left. 

If there are no intersections, it means that no values 
of C,; and Cp at the chosen speed and acceleration can 
be found. 

If a point of intersection of a circle of the form of Eq. 
(8) with the polar curve is found, three cases may be 
represented: 

(1) When the value of Cp in the point of intersection 
of circle and polar curve is greater than the value of 
Cp in the center of the circle, being (An(hp.)/Vq¢f) — 
(W/gqf)(dV/dt), the requirements of Eq. (7) have not 
been met with; 
flight. 

(2) When the value of Cp in the point of intersection 
is equal to the value mentioned in Case (1), Eq. (7) 


the speed is too high for continuous 
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is satisfied only by putting L = Wor C, = Cj,. There 
can be no more acceleration because the plane is flying 
horizontally with the throttle entirely open. 

(3) If the value of Cp in the point of intersection is 
smaller than the value mentioned in Case (1), flying is 
possible and an excess of engine power is available for 
climbing and acceleration, since, in this case, the co- 
ordinates of the point of intersection of circle and polar 
curve satisfy Eq. (7). 

If the drag coefficient Cp is known, it is a simple 
matter to determine the climbing speed C = V sin @ in 
the usual manner. 


? 


Sin @ can also be read from Fig. 3. 
An(hp.) 


; aa (* (a) = C+ C 
Vaf geqf/\ dt, 


In this equation (An (hp.)/Vqf) — (W/gqf) (dV/dt) 
is regarded as the total available drag coefficient for 


(13) 


flying and climbing and, since Cp is the drag coefficient 
of the plane corresponding with the speed V and the 
climbing angle @, Cc can be regarded as the drag coeffi- 
cient available for climbing, so that for the rate of climb 


CW = Coa fV or C = V(Co/Cu) 
and 
sin @ = Co Cra 


The angle @ can be immediately read from the figure. 


DISCUSSION 


(1) Considering Eq. (9), it is clear that the place of 
the center of a circle is dependent on two terms, of 
which Ay(hp.)/Vgf in a chosen speed is known but the 
second term is not. It can be assumed that the best 
rate of climb occurs when D/C, is constant and when 
q is constant, so that d\’/di can be found as a function 
of pand?. This, however, does not give an easy solu- 
tion because, before kno xing C as a function of ¢, it is 
difficult to express p as a function of ¢. 

However, the approximate rate of climb can be de- 
termined, assuming that there will be no acceleration, 
in which case V can be found as a function of ¢, which 
is sufficiently satisfactory for the construction of the 
real curve of the rate of climb. 

(2) In the preceding consideration a chosen speed 
is always mentioned. This chosen speed is a function 
of the Cp — C, curve and the propulsive efficiency. 
It is evitent, however, that the best rate of climb will 
only occur when flying with the same airplane but with 
an assumed weight Z at a speed that does not much 
differ from ‘the speed with a minimum effective pro- 
peller horsepower, in which case Cp/C,* * is a minimum 

If nis known or sufficiently valued, the speed in ques- 
tion, with minimum horsepower, can be found ap- 
proximately from the equations 


'/, Cy V*f = G cos @ and 


tan 6 = [An(hp.)/Vq@fC,] — (Cp/Cz’ 
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Fic. 2. Construction of Ce—no acceleration assumed. 


and the values of C, and Cp, for which Cp/C,¥? is a 
minimum. The speed found in this way has to be cor- 
rected with respect to the course of the propulsive 


efficiency. 


PROCEDURE 


(1) Determine the best speed.for climb. 
(2) Plot in the C,;—-Cp curve for the whole airplane 


the point: 
_ An(hp.) ) 


Vof 

(3) Determine the value of Cz,,, according to the 
speed V, the air mass density p and the weight of the 
airplane W. 

(4) Describe a circle with Cz), as the radius and the 
point determined in item (2) as a center. 

(5) Determine the point of intersection of the circle 
and the C;-—Cp curve. 

(6) Calculate from the figure the value of C¢ ac- 


C, = 0 (see Fig. 2) 
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Fic. 3. Construction of Ce—acceleration assumed. 


cording to Eq. (13) assuming that there is no accelera- 
tion. 

(7) To find the climbing speed C, multiply Cc by V 
and divide by Cz,. 

(8) Perform this at the various intervals of heights. 
Correlate in each interval the rate of climb according to 


the acceleration. 
(9) Determine, with the help of the values found, 


the value of (W/gqf) (dV /dt). 
(10) Plot in the C,-Cp curve the point 


= Arte) _ (WY), 6, oo eerie. 
Cp Vaf — Cr = 0 (see Fig. 3) 


(11) Continue with the same process as described in 
items (4) through (7), with the point determined in 
item (10) as the center of the circle. Note the value 


— An(hp.) 
Vof 
(Continued on page ea 





Space Limitations and Optimum Conditions 
in Aircraft Spring Design 
RICHARD H. CARTER* 
United Aircraft Corporation 


SUMMARY 


The effect of space limitations on spring performance is analyzed 
and leads to optimum design formulas for the three most com- 
mon types of aircraft springs—i.e., tension, compression, and 
torsion springs. Using the formulas thus derived, for a given 
allowable stress the relationship between the space limitations, 
the maximum load obtainable, and the corresponding wire size 
can be found immediately. In addition, the process of spring 
design becomes direct and eliminates trial-and-error methods. 


GENERAL 


Ir PRACTICAL AIRCRAFT spring design, the designer is 
asked to specify a spring that will fulfill certain re- 
quirements. These requirements may vary somewhat 
but are usually ramifications of the general problem 
that the spring must do a certain job yet conform to the 
available space. The procedure usually employed is 
“trialanderror.’’ The designer uses the various formulas 
and charts and by intelligent guessing tries to find a 
combination of the variables which will meet the re- 
quirements. But often, when the spring requirement 
exceeds, or is just within, design possibilities, this is a 
tedious process. 

To improve this procedure, certain optimum condi- 
tions are developed in this paper. Given a limited 
space, the use of these conditions affords direct calcu- 
lation of maximum performance. Thus, if the require- 
ments are too stringent for proper spring design, this 
may be quickly detected. And as the most efficient 
spring design is given by these formulas, any surplus 
weight or space caused by nonoptimum design is 
eliminated. 

Besides the applications to efficient spring calcula- 
tion, the airplane designer who wishes to incorporate a 
spring in a unit can immediately calculate the minimum 
space for proper spring operation without actually hav- 
ing a spring designed. This should eliminate the com- 
plications that arise later when it is found that inade- 
quate space has been left for the spring. 

To develop these optimum conditions, the maximum 
load as determined by the allowable stress formula is 
equated to the load as determined by the formula in 
terms of the physical characteristics of the problem (i.e., 
the outside diameter, maximum deflection, etc.). Vari- 
ous manipulations are then performed to obtain the re- 
sults in the desired forms. The accompanying align- 
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ment charts are constructed for stresses normally en- 
countered in aircraft spring design and reduce the cal- 
culations to a minimum. 

Although the optimum design conditions for the most 
common general cases are presented in the main analy- 
sis, particular problems with various modifications may 
be solved by logical extensions of the formulas. Im- 
portant examples at the end of the paper indicate this 
procedure; the last example is particularly important 
since it indicates the general approach to certain prob- 
lems in spring design where the weight must be a 
minimum. 


EXTENSION SPRING ANALYSIS 
The formulas usually used in spring design are 


Pmar, = 47d°S/8K,D 


and 
Prog = [Gd*/8nD* 


where 


c — 
5, = Wilt elec = 1 + o.o1()( 22 — 414) 
D/\ D 


S = stress in wire 

G = torsional modulus of elasticity 

D = mean diameter of coils (D, = outside diameter) 
h solid heighth of coils 

d wire diameter 

f deflection 

nm = number of coils 


Equating Pmaz, and P,,,, introduces optimum condi- 
tions on the variables; solving for d gives 


d? = SxnD*/K,fG (1) 


To eliminate and introduce h, the optimum relation 
that n = h/d is used; Eq. (1) then becomes 
d = DV Sah/KifG (2) 


K, is a function of D and d but is too complicated to 
handle conveniently. Hence, an approximation, K,’ = 
[5D/(5D — 3d)]*, is used. (K,’ is accurate within 1 
per cent in the range 4 < (D/d) < ~.) 

Substituting K,’ for K, in Eq. (2) gives 


d = 5D/(3 + 5VGf/Sah) (3) 
Substituting (D, — d) for D in Eq. (3) gives 
d = D,/(1.6 + WGf/Sth) (4) 
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(This is the optimum design wire size for a given f, h, 
and D,.) 

d is now eliminated from the given formula (Pima: = 
aSd*/8K,D) and Eq. (4), giving 


ro (MMe as) © 
maz, ~~ ‘4 ra * 
8 /\KiD/\1.6 + VGf/xSh 


But KD is approximated within 2 per cent in the 
complete range from 2 < (D/d) < ~ by D,/[1—(d/2D,)]. 


ra (2 casvarewil 8) © 
mez, ~ 8 (1.6 + V Gf/Sh)? 2D, ), 


But 1 — (d/2D,), with d = D,/(1.6 + WGf/xSh), is 


equal to 
(1.1 + WGf/wSh)/(1.6 + © Gf/rSh) (7) 

By substituting Eq. (7) in Eq. (6), 

wSD (11 + VGf/rSh) 


ee (8) 
8(1.6 + WGf/nSh)4 


P nas. = 
(This is the maximum load possible for a given f, 4, and 
D,.) 

If the designer is given the spring rate—i.e., P at |; 
and Pyar at ,—and the minimum space for proper 
operation is desired, it is desirable to have explicit ex- 
pressions for D,, f, h, and d in terms of the given values. 
To retain reasonably simple expressions for these 
values, an approximation to Eq. (6) is made by assum- 
ing 1 — (d/2D,) equal to 0.941 (this assumption is ac- 
curate within 3 per cent in the range 4 < (D/d) <= 
and exact for D/d = 7'/2). 

With this assumption, 

/8.5)(1.6 + V Gf/Sh)-* 


Fs Zz. — (rSD,?, (8a ) 


Solving Eq. (8a) for f/h gives the maximum value for 
f/h for any D,. 
aSD,? 


aS SD,2 \"/* 2 
l/h = — ——} — 16 9 
ffl G (4 =) | (9) 





| ae | 


f =(h -—h)/0 — (P/ 
Or solving Eq. (9) for D, gives the minimum value of 
D, for any f/h. 


8. OP maz. ) ‘| ; 
Dew | es 1.6 
Ee (y¥% sila 


And from the given equation Png, = md°S/8K,D, 


d = (See a (1 1) 
7S 


(10) 


But 
(K,D)'* = D,' (Kil — d/D,)\'” 


or approximately, 


(KD) 


= 0.975D,[1 + 0.41(d/D,)] 


TICAL SCIE 
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Substituting this in Eq. (11) and solving for d gives 
the approximate expression 


1 i 4 iE 4. . (?: Fass.) | (11a) 
a D.\ Ss 


(The accuracy is within 1 per cent in the range 2 < 
(D/d) < @.) 
Or approximately, 


d = (8.5D,P maz /'®S) (11b) 


(exact for D/d = 71/2). 


COMPRESSION SPRING ANALYSIS 


The same analysis as for the extension springs is ap- 
plicable, except that the two inactive end coils change 
the relation m = h/d. Usually the extra length is con- 
sidered as 1!/.d, but various considerations have shown 
it more practicable to use 2d in this analysis. 

Calling h — 2d = h’, from Eq. (4) 


= D,/(1.6 + VGf/rSh’) 
But 


fr _ oii 1 oy 
Wi” Vi F =e cars | 
or very nearly 


V ish’ = [1 + (d/h) V1 /h 


Therefore, 
D, 
1 = sanyo 12) 
1.6 + V Gf/rSh[1 + (d/h)| 

It would be possible to solve this quadratic for d, 
but a simpler expression and one analogous to that 
in the extension spring analysis is found by putting in 
the approximate expression for d in the right-hand mem- 
ber. Eq. (2) is used with nm = h/d, D = D,, anda 
representative value of K = 1.41; these approxima- 
tions will introduce no practical error for d in the right- 
hand member. Thus, 





D 
0 


; 13 
1.6 + VGf/xSh + (D,/1.2h) aad 


d= 
And by the same procedure as in the extension spring 
analysis, 
aSD,” {. 1 + Vv ‘Gf/mSh + - (D,/ 1. 2h)] 


—— (14) 


P se. = 
8[1.6 + VGf/xSh + (D,/1. 2h)|§ 


From the derivation, it is easily seen that Eqs. (11), 
(lla), and (11b) are applicable to compression springs 
without modification. The other expressions for D, 
and f/h may be derived, but it is simpler to use the ex- 
tension spring formulas with the correction that 
V Gf/xSh i in the extension spring formula is equivalent 
to V Gf/2 rSh + (D,/1.2h) for compression springs. 















AIRCRAFT 


TORSION SPRING ANALYSIS 


The ordinary formulas are 
aM - ( 6 \— Ed* -) 
ron peo 25Dn 


= inside diameter of 


3 


WU oe. 
4" mar, ~ 10.2K 


mean diameter of coils (D 
coils) 
= number of coils 
= Wahl’s factor for torsion springs 
wire diameter 
angle of deflection 
stress in wire 
= Young’s modulus 


The correct wire size for any maximum moment is 
given immediately by solving the first of the above 
Thus, 


d= (10.2K M maz,/S)' 7 


equations for d. 
(15) 


therefore 


But K’’’ is approximately [1 + (d/5D,)], 


10.2M maz d ) 
j« (= ; 
aaa taae ar 


or solving for d gives the approximate expression (ac- 
curate within 2 per cent from 2 < (D;,/d) <=) 


wa | 1 | | 
d = 1+ “aes 15: 
( 5 BDA. S oe 


Or approximately, 


2.25( M maz. /S) ’ (15b) 


(exact for D;/d = 6). 


To find the optimum wire size as a function of the 
available space, Wn; is equated to J/, in the given 
formulas. Solving for d gives d = 4,050DSn/10.2K9. 

To find the other optimum relationships, the optimum 
condition that n = h/d is again used. 

But D/K = (D; + d)/K. Since K is approximately 
0.965 [(D; + d)/D;|, D/K = D,/0.965. Making these 
two substitutions gives 


4,200S / Dd 
d? = 102E ( *) or d = 20. 3y(5 ve “ (16) 


(The error in this assumption is less than 2 per cent 
in the interval 3 < (D;/d) < ~.) 

Eliminating d between Eqs. (15) and (16) gives the 
equation relating the maximum moment available with 
a given space and angle of deflection. Thus, 


fo is (42) ( S =) (2#)" 
a 10.2 6 


1/K is approximately 1 — */;(d/D,) and, by substi- 
tuting Eq. (16) for d, Eq. (17) becomes 
(4) {o 7a) 


5 3 }/2 
Muer = 819) = (74)"] [1 - 29 
, ms 8 


(17) 
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(This approximation is accurate within 1 per cent in 
the range 3 < (D;/d) < @.) 
Or approximately, 


S\*/Dih\' 
we a 7 50( ) ( ) 
E 8 


SPECIAL PROBLEMS SOLVABLE BY  UR- 
THER ANALYSIS 


EXAMPLES OF 


Most spring problems are easily transformed into a 
form easily solvable by the formulas developed in the 
general theory, and, with a little experience in their use, 
the designer can quickly determine the most efficient 
spring or talculate the space versus performance rela- 
tion. 

However, there are some problems of a more compli- 
cated nature which are particularly important to the 
aircraft spring designer. These are readily solved by 
extensions of the general principle of optimum design 
as illustrated by the following examples. 

These three examples are for extension spring; how- 
ever, the methods of analysis are quite similar for the 
other two types and give analogous results. 


Example 1 
Given: 
Allowable outside diameter (D,). 
Load P» at length /s. 
Length J,(l; > 2). 
Length required for end connections (e). 
To find: 
Lowest load (P; = Paz.) at hh. 
Corresponding wire diameter (d). 
From the ordinary formulas, 
Pose = 90°S/8K sD 
fGd* 


} ws ss a 
ShD* 


fGd® 
8D*(l, — f — e) 
Eq. (2) introduces two optimum conditions: (1) h = 
nd; (2)ht+e=1]-—f. 
The third relation between P,,,, and P, is obtained 
from the spring-rate expression : 
P nas. f _ P, (f + ls ae lh) (3) 


Eliminating P,,,2, and f between Eqs. (1), (2), and (3) 


gives the equation for the wire size: 
d* hL—e 
8P:.D,—d) 7S | Gd*(l, — 1s) 
(ly —e 
oo (D, — d)* 


-=Q (4) 


Although Eq. (4) is a fifth degree equation in d, it 
is easily solvable by an approximate method, since a 
reasonably close answer for d may be estimated from 
the conditions of the problem and the use of the opti- 
mum design charts. The equation contains no ap- 
proximations and, thus, an answer to any degree of ac- 
curacy may be obtained. The author has found it ex- 
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OPTIMUM DESIGN OF EXTENSION OR COMPRESSION SPRINGS 
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pedient to set up forms for computing residuals for the 
selected values of d and then roughly plot these residuals 
against d; where the residual is zero gives a close ap- 
proximation for d. Thus, the correct wire diameter is 
easily found in about 15 min., and the designer knows 
that the answer is the best possible solution to the 
problem. 

An explicit expression for P,gz, would be complicated, 
so the quickest way to find Pg, is by substitution of 
the optimum wire diameter in Eq. (1). 


Example 2 
The result of the following analysis will often be of 
use as an optimum design condition. 
Given: 
Allowable outside diameter (D,). 
Length /». 
Length /, (J; > 2). 
Length required for end connections (e). 
To find: 
The wire diameter (d) which will give the maximum 
load, P., at /, without overstressing the spring 
at /,, and the corresponding total extension (f). 


Solving for P, in Eq. (4) gives 


l 
S(1,; —e) 


1S d3 : dé 
— (le — e) -— G(, — L)——_—— (4a) 
Ky, D, —d (D, — d)8 


The maximum value of P, is found by setting the 
derivative of P2 with respect to d equal to zero. 
Wahl’s factor (K,) must be taken into account since 


it is a function of d. Using 


3( d bee : 7 
1+ 0.015 ( 
D,-—-d D,—d 


OK, 1.4D 


K, = 
gives 


Od (D, — 2d)(D, — d) 
to a good approximation. 


Using this value of OK,/Od in the expression for 
OP,/O0d = 0 gives 


0.64S(l. — e) ae Ky “ 
G(l; — 12) ((D,/d) — 1}? 
| (D,/d) — 0.4 (5) 
0.465 


)P. d — 0.667 — — \ 
K,[1 — (2d/D,)| 

The values of (J; — /,)/(l2 — e) (see Fig. 1) are readily 
plotted against D,/d giving a convenient method of 
calculation. . 

A good approximation for the value of f to give the 
maximum P, without calculating the wire diameter 
may be found by eliminating d and D from Eqs. (1), 
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(2), and (3) giving P. in terms of f. Assuming K, a 


constant and setting 0P2/df = 
iu (1, — 5) (1; — ¢7 (6) 
(1; — 0.4/4, — 0.6e) 


0 gives 


the value of f obtained being somewhat high. 


Example 3—Optimum Design with Respect to Weight 

In Example 2, the wire size, solving the given prob- 
lem with the specified space limitations, is found by 
Eq. (4). 
limitations and not the space limitations are critical. 


Often in aircraft spring design, the weight 


In many cases there is an optimum design with respect 
to the weight, as is illustrated by the following exarr ple. 
Given: 

Load P, at length /. 

Length /,(/; > ds). 

Length required for end connections (e). 

To find: 

Lowest load (P; = Pna;.) at 1; consistent with the 
corresponding wire diameter (d) and coil diam 
eter (D) combination giving the lowest weight. 

From Eq. (4), (D, — d = D) 


r |= — e) G(l, — 1,)d* 
da’ _ 


—8P,D(h—e) (7) 
K, D? 


But the weight of a spring (neglecting the inactive 
portion) is proportional to the cross section of the wire 


multiplied by the length. Or, 


W = (constant)d*7D = (constant)dhD(n = h/d) 


Therefore, the quantity (dD) is to be minimized 
(since / is a constant) where d and D are restricted by 
Eq. (7). 








AIRCRAFT 


Using the Lagrangian multiplier (A) method to ob- 
tain this minimum value of dD and using 


OK,/dd = 1.4/(D, — d) 


0K,/0D = 1.4d*/D(D — d) 


as approximations for the derivatives of K,, the two 


equations for \ become 
3mS(I2 — €) (1 __ 0.467d ) mi 
K,(D — d) 
5d°G (1; =) 7 
D2 





r 


D + dd? | 


1 


ie a! aad ais 
t+ mn} -22 mS(ly — e)d* , Gh 1d 2 
K,*(D — d)D Ds 


4P,(]; — o)| 


Eliminating \ and calling A = (mS/K,)(/, — e), 
B = 2G(l, — l), and C = 16P2(l, — e) gives the equa- 
tion 


7B(d/D)* — C(D/d*) = 6AR (8) 


(Continued from page 50 


As may be seen from the foregoing, the method de- 
scribed is extremely simple but, assuming that the 
C,-Cp curve for the whole airplane is drawn accurately, 
more precise than the usual methods. With most of 
these methods, continuation of the calculation of the 
rate of climb curves at speeds under the minimum hori- 
zontal speed of the airplane is not possible. And, as 
the calculations show, it is possible to fly under the so- 
called minimum speed when climbing. The way in 
which the acceleration is taken into account is only a 


SPRING 


DESIGN 


where R is nearly equal to 1 — 0.9(d/D). 
Rewriting Eq. (7) gives 
B(d/D)? + C(D/d*) = 2A (9) 


Eqs. (8) and (9) are linear in (d/D)* and D/d'; 
solving for these quantities gives 


(d/D)* = (A/B)[{1 — (0.7d/D)}| 


D/d* = (A/C){1 + (0.7d/D)] (80) 


To find the value of D/d, a plot is made (see Fig. 1) 
of D,/d against (1; — k)/(l — e). With this value, 
d is found from the combination of Eqs. (10); i.e., 
from 


d? = (C/A)(D/d)/[1 + (0.7d/D)] (11) 


THE ALIGNMENT CHARTS 


These charts are primarily intended for use as an aid 
in estimates and as a check on the formula calculations. 
However, as all but the most exact spring design prob 
lems do not require a design accuracy within 5 per cent, 
the charts alone are usually sufficient for calculation. 
In the author’s own experience, these have been a great 
aid in nearly every spring design problem encountered 


more exact repetition of the work already done in 
determining the rate of climb in a first approximation, 
which is already more accurate than that followed in 
most of the usual methods. 
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Dynamics of Constant-Speed Propellers 
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SUMMARY equations of motion for the controlled propeller will be 





A theoretic investigation of the performance of constant-speed developed. Let 
propeller controls, operating with various types of anticipatory 

















components and several classes of lag, is undertaken over the A = propeller speed (r.p.s.) 
performance range of the airplane. The effect of aerodynamic n = change in propeller speed from the de- 
parameters of the propeller on the performance of the overall sired speed (r.p.s.) 
system is established. The advantage of nonlinear control 8 = blade angle (degrees) 
characteristics is noted and demonstrated by calculation of . \ 
; ane : Bi = change in blade angle (degrees) 
control response for a typical combination of propeller, engine, : baie, iia h.) 
control, and airplane. A number of charts are presented in J = forward speed of the airplane (m.p. 1. | 
dimensionless form, giving the propeller derivatives involved in v = change in forward speed of the airplane 
control design. Charts showing stability requirements for a (m.p.h.) 
control system subject to both first- and second-order lag are QO, 0., Q, = torque, engine torque, propeller torque 
supplied, as are charts for stability determination of control (ft.1b.) 
systems with constant time lag. engi : 
AQ = change in torque (ft.lb.) 
INTRODUCTION I = inertia of the engine-propeller combina- 
tion, referred to the propeller (Ib.-ft.- 
I IS DESIRABLE to maintain constant propeller and sec.?) 
engine speed for airplanes under changing conditions 





A convention of signs will be chosen by defining a 
positive torque as one that tends to accelerate the 
propeller. The acceleration of the propeller is then 
proportional to the algebraic sum of engine and aero- 






of flight.!. This regulation is obtained in conventional 
designs by varying propeller blade angle according to 
some function of the speed error. The control consists 
of an element capable of measuring the speed error and 
its integral or derivatives and a servomechanism capable 
of moving the propeller blade in accordance with the 
indication of the error sensitive element. Design de- 
tails of successful controls are available in the litera- 
ture.” * 







dynamic torques. 
Engine torque is a function of several variables— 






among them, engine speed. Increase of engine torque 
with speed is greatest on engines with high blower 


























































The problem of control design is somewhat com- ” es 
plicated by the variation of the aerodynamic parame- eet 
ters of the propeller over the operating range of the # | 
control. Since it is customary to keep the sensitivity - ‘ pT LYwriGHT R-1820-63 
of the control invariant with airplane velocity, altitude, | gia is 
and engine torque, while change in propeller torque es ae | | | || 
with blade angle and propeller speed varies with for- 0, 
ward velocity, density, and power coefficient, it is 22 | 
possible for a control that is perfectly satisfactory for aoe 
one flight condition to hunt undesirably at another. feo = 0 —— oe } : 
Although the propeller and control constitute a non- | ve eee a | | 
linear system with coupling ‘factors, damping and -" — T we Woouece wase" | : 
spring coefficients functions of the independent vari- | ite! “aoa” 
ables, it is possible to assume constancy of the principal 2" | TY _ CRUISING a 
parameters by restricting control motion to small ©° | | 
displacements. Response characteristics can then be . 4 “S00, | cals oot Tn, 
determined by separately examining points distributed p ERR TREE Ce,” * = Se proves’ — aa | ’ 
over the region of use of the control. }. Smee sees | 
| i, ee | 
10 0 FT. = = Ss _ . | 4 
PROPELLER EQuaTIONS = | | oe | 
; {| [| _ /} cruising ' ' 
Assuming that regulation of propeller speed is to be * oon, | eoooort.y fT LYCOMING R-680 | 
obtained by control of propeller blade angle, the P ] = et ON EE seh 
‘ 1400 1600 1000 2000 2200 2400 


Received August 18, 1942. 
* Antiaircraft Artillery Board. Fic. 1. Torque variation with speed of typical engines. 





DYNAMICS OF 


ratios; unsupercharged engines develop substantially 
constant torque with changing speed over their normal 
range of operation. Torque variation for a number of 
commercially available engines is shown in Fig. l, 
where an increase of about 0.80 ft.lb. per r.p.m. is 
noted for a Double Wasp S1A4-G. Data on high- 
powered military engines where torque gradient with 
speed is expected to be greatest are, of course, unavail- 
able for publication. Change in engine torque can be 
expressed, 


AQ. = (02,/0N)n + AQo (1) 


where AQ, includes changes in engine torque from 
other sources than engine speed, such as throttle 


setting. 
For the propeller 


Q, = f(N, V, B) (2) 
and for small changes in the independent variabies 
00, 00, 00, 
AQ, = = — —v 3 
Q, aN 2 8 Bi + Pa) r ( ) 
Acceleration of the propeller is given by 


AQ. + AQ, (4) 


In = 


oN 


) ), » a, : 
= (22 + Cn a 7 Bi + - ty + AQ (5) 


op. 3) 

The dimensions of Eq. (5) are correct if m is expressed 
in radians per second. However, it will be more con- 
venient in numerical work to consider m in revolutions 
per second, in which case inertia J is replaced by 


2rI. 
The following shortened notation is then adopted, 


Qyv = 1/2nI[(0Q,/0N) + (00,/0N)] (6) 

Qs = 1/2rI(0Q,/08) (7) 

Qy = 1/2xI(0Q,/0V) (8) 

Qo = 1/2nI(AQo) (9) 

and the equation of motion of the propeller becomes 

n — Qyn — Qghi = Qo + Ovo 

A “time constant” for the propeller is defined by 

assuming the blade angle to be fixed and noting the 

response of the propeller to a torque Qo suddenly ap- 

plied when the propeller is turning at some constant 

speed N. The change in propeller speed after the 
initial transient is 


(10) 


1, = —Q)/Qw (11) 


and the propeller reaches this speed change in an ex- 
ponential approach 
1 — ¢-¥% 


(12) 


n/ny, = 


where 


T, = —1/Qy = propeller time constant (13) 


CONSTANT-SPEED PROPELLERS 




















Fic. 2. Aerodynamic torque gradients. 


AERODYNAMIC TORQUE DERIVATIVES 


The aerodynamic torque derivatives Qy and Q, can 
be calculated for any particular propeller by plotting 
torque against propeller speed and blade angle in the 
neighborhood of the point under consideration and 
noting the respective slopes of the curves. This has 
been done for a particular case, and the results are 
plotted in Fig. 2. The engine is a Wright R-1820 G-3, 
geared 16:11, delivering 840 b.hp. at 2,100 r.p.m. and, 
8,700 ft. Propeller diameter is 12.0 ft.; inertia of the 
propeller-engine combination referred to the propeller 
shaft is 1,300 Ib.-ft.? or 40.4 Ib.-ft.-sec.*. Air speed is 
250 m.p.h. Then by Eq. (6), the propeller’s con- 
tribution to Qy is 


Qnp = + ( )(—9.0)(60)(16 11) = —3.09 per sec. 


Contribution of the engine torque gradient is 


1 acs 
On = + (=. ).o.s0y(60)(16, 11) =+0.28 per sec. 


From Eq. (7) 
1 , 
Og = + (+) (—600) = —2.36 per deg.-sec.? 
ot ° 


A more general approach to the torque derivatives is 
through the conventional propeller coefficients. The 
propeller torque can be expressed in terms of a power 
coefficient c, as 
(14) 
where p is air density, D is propeller diameter, and c, is 
available from any of the standard references on pro- 
peller data. Differentiating Q, with respect to pro- 
peller speed and blade angle and making use of the 
following notation 


] Oc, 
on = 3, (2 - 735) 


pe = (1/2) (Oc,/08) 


Q, = cpN*D8/2x 


(15) 


(16) 
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y= (1/2z1) (pszD*) (17) 


where ps; is standard air density at sea level and 
J = V/ND, it is found that 


Ov = qvnNo (18) 
Os = ganN2o (19) 


Here g, and gg, are dimensionless parameters de- 
pending only on the aerodynamic characteristics of 
the propeller. 7 is a dimensionless parameter depend- 
ing only on propeller diameter and effective propeller 
inertia. o is the ratio of actual air density to standard 
air density at sea level. Fig. 3 shows n computed for a 
number of propellers of various diameters from data 
found in propeller literature‘ or supplied by the indus- 
try. Engine inertia has not been included in the 
computations for y in Fig. 3, and a proper increment 
should be added to account for engine inertia in any 
particular computation where it is known. 
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Fic. 3. Propeller inertia factors. 


Computation of cy and gg has been carried out for two 
typical propellers with blades of Clark Y section, a two- 
bladed propeller’ and a three-bladed propeller.® Fig. 4 
shows the variation of cy with J and blade angle. Since 
propeller torque increases with propeller speed, oppos- 
ing acceleration, cy is negative. It is noted that the 
absolute magnitude of cy increases rapidly with blade 
angle and more slowly with J. gg is shown in Fig. 5. 
dg also is negative, increasing in absolute magnitude 
with J. Variation of gg with blade angle, however, is 
neither large nor distinct. 


LAWS OF CONTROL 


In order to secure control of propeller speed, the blade 
angle is varied according to some function of speed 
error. Change in blade angle in turn causes change 
in propeller torque. If the processes involved between 
measurement of the speed error and movement of the 
propeller blade can be expressed in the form of linear 
differential equations with constant coefficients, the 
control is considered to be ‘“‘linear.”” Actually, an 
exactly linear control is a rare exception. However, 
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Fic. 4. Variation of gy with J and blade angle. 
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Fic. 5. Variation of gg with J and blade angle. 


the assumption of linearity can usually be made, sub- 
ject to whatever restriction to small movements of the 
control may be required, in order to obtain with readily 
available mathematical tools results closely if not 
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exactly, definitive of the actual behavior of the con- 
trol. 

Probably the simplest law of propeller control would 
be that making the change in blade angle from equilib- 
rium proportional to the speed error. This type of 
control has been treated at some length in an earlier 
paper by the writer.’ It has the disadvantage that 
with the control adjusted to give the desired -propeller 
speed for some particular engine torque and air speed 
any departure from these conditions will cause an 
error in propeller speed which will be present until 
torque and air speed return to their original values. 

By causing the change in blade angle to be propor- 
tional to the time integral of the error in propeller 
speed, the control is given the ability to wipe out steady 
state errors arising from constant departures of engine 
torque and air speed from particular reference values. 
It is frequently further desirable to make the change 
in blade angle proportional to the derivatives of speed 
error. 

Then assuming temporarily that there is no lag in the 
control, it is possible to write the most general ex- 
pression defining response of a linear control as 


B=k fndt+kn+kn+... (20) 
where 
ky = 2B:/dn 
ke = O8;/O0n 
ks = CBi/Cn 
The physical nature of the integral control will be 
made more evident by differentiating Eq. (20) once 
and setting k, and higher order components equal to 
zero. Then 


dB/dt = kn (24) 


from which it is seen that the integral component makes 
the rate of change of propeller blade angle proportional 
to speed error. The Hamilton constant speed control? 
is of this type. 


CONTROLLED PROPELLER 


Upon combination of Eqs. (10) and (20), the follow- 
ing expression results. 

i(1 — ksQg) — (Qn + keQg) — kiQgn = Qo + 8Qy (25) 

The following observations can now be made regard- 
ing speed error without the necessity for solving Eq. 
(25). 

(1) For constant change in torque Qo or constant 
change in air speed vw, there is no steady state error in 
propeller speed. 

(2) For constant rate of change of torque Qo, there is 
a steady state error in propeller speed given by 

Nss = —Oo/kiQe (26) 

(3) For constant acceleration dv/dt along the flight 
path, there is a steady state error in propeller speed 
given by 


CONSTANT-SPEED 
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— (dv/dt)Qy/kiQe (27) 


Ns, = 


It is convenient to replace —Q,/Qg¢ in Eq. (27) by its 
equivalent from Eqs. (7) and (8) 


—Qy/Qs = (08/0J)/(ND) 


(28) 


whence 


— (dv/dt)(08/0J)/(NDk;) (29) 


Ny; = 
Parameter 08/0/ depends only on the aerodynamic 
characteristics of the propeller. It will be mentioned 
again later in connection with maximum required 
blade rate, at which time charts for its determination 
will be provided. 

(4) Damping for the motion represented by Eq. (25) 
consists of aerodynamic damping from the propeller 
(—Qy) and ‘‘control damping,” supplied by the control 
component that makes change in blade angle propor- 
tional to speed error. As indicated by the numerical 
computation of Qy, an increase in engine torque with 
engine speed makes the numerical value of Qy less 
negative, and consequently the damping of the system, 
which is — Qy is reduced. 

(5) Effective inertia of the propeller is reduced by 
the component of the control which makes blade angle 
proportional to the first derivative of the speed error. 
By making k3Q, = 1.00 the effective inertia of the 
propeller is made zero. Actually, lag in a physical 
control would prevent this ideal from being attained, 
and the control would become unstable, tending to 
hunt as k;Qg, approached 1.00. 


DAMPING RATIO 


The simplest type of integral control is that which has 
only a k; component and depends for its damping on 
the aerodynamic damping of the propeller. Defining 
damping ratio {) as the ratio of actual damping to the 
minimum amount required to produce aperiodic motion 
of the system after a disturbance, it follows from Eq. 

25) that for the simple control 


fo? = —Qy?/4kiV, (30) 


Converting Eq. (30) by substituting gy and gg from 
Eqs. (18) and (19), 


fo? = —(qn?/4q@)(no/hi) 


Damping ratio is then decreased by: 

(1) Smallo. As altitude increases, the control tends 
to become more oscillatory. 

(2) Small. Increasing propeller inertia for a given 
diameter or decreasing propeller diameter for a given 
inertia will cause the control to become more oscillatory. 

(3) Large k;. A sensitive control, causing large 
change in blade rate for small speed error, tends to 
cause oscillatory motion. 

(4) Small value of vy, where 


v= V —qn?/4q0 
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This ratio involves only general aerodynamic char- 
acteristics of the propeller and has been computed for 
two propellers. Fig. 6 shows v vs. J for various con- 
stant values of O0C,/OJ. It was considered desirable 
to choose 0C,/O0J as a parameter, since, for constant 
propeller speed and constant engine torque, 0C,/0J 
is constant as the airplane changes its forward velocity. 
This is evident from the expression for the power 
coefficient C, 


C, = V¥/p/PN (33) 


where P is horsepower. 
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Fic. 6. Variation of v with J and blade angle. 


Change in airplane speed for constant power moves 
the airplane along a curve of constant 0C,/O0/. Change 
in power for a constant air speed moves the airplane 
vertically along a line of constant J. The following 
further conclusions regarding damping of the motion 
of the controlled propeller can then be drawn with the 
aid of Fig. 6. i 

(5) Damping decreases with decreased forward 
velocity for given engine power and r.p.m. 

(6) Damping increases with increased power at a 
given air speed and r.p.m. for constant engine torque 
gradient 00,/ON. 

Since Fig. 6 indicates that aerodynamic damping is 
least at zero velocity of the airplane or J = 0, this 
region is further investigated in Fig. 7 for propellers 
of two, three, and four blades.’ Damping increases 
with number of blades and, as might be expected, with 
blade angle. 

The curves of v in Fig. 7, if continued to smaller 
angles 8, would reach a minimum and again rise, in- 
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Fic. 8. Variation of 08/0J with J and blade angle. 


creasing to infinity since —gy is always positive; but 
—qg passes through zero at the blade angle of minimum 
torque when a change of blade angle in either positive 
or negative direction causes an increase in torque. 
In this region the assumption of linearity of the system 
for small movements no longer holds and, indeed, 
control in this region is not possible by the system dis- 
cussed in this paper. 


MAXIMUM BLADE RATE 


The maximum airplane acceleration with which a 
control must cope is one of the factors determining the 
maximum rate of change of blade angle which the 
control must provide. Disregarding transient re- 
quirements, it follows from Eq. (10), when propeller 
speed error is held negligible, that 


Bmaz- = (dv/dt)mar.(08/0J)(1/ND) 3 
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Again 08/0J depends only on the aerodynamic char- 
acteristics of the propeller and is plotted in Fig. 8 for 
two propellers. It is noted that for a given propeller 
speed, engine power, and airplane acceleration the 
maximum blade rate is required as the airplane passes 
through a region at about J = 1.30. 


CONTROL RESPONSE 
The following notation is now introduced 
5 = k2Qe/Qn 
vy = —k3Qg = (1/D(0Q/dn) 


Here 6 is the ratio of artificially introduced damping 
to aerodynamic damping already present in the con- 
trolled system, and y is the ratio of change in control- 
ling torque for a unit change in propeller acceleration 
to change in system torque (inertia) for a unit change 
in propeller acceleration. 

Eq. (25) is now 


(35) 
(36) 


n 
P23 


1 1 
(1+) +2048) +3 ( 1) = +90, (37) 
D Pp t80 


The damping ratio of the system with derivative 
controls is 


42 = fo%(1 + 8)?/(1 + 7) (38) 


and the solution of Eq. (38) for the error in speed, 
following a small sudden change in torque AQo, is 


—§ 


a 
” = Ae ( +7) sin w (t/T>) 


(39) 


Ny, 
where 


A = [1/(1 + 8)] (26./V 1 — §,?) (40) 


w = [(1+8)/1+ YI(V1 — 2/28.) (41) 


For the case of aperiodic motion of the system, when 
¢, = 1.00, 

n 7 ( 1+é 
= = [1/1 + MU/Ty et 

It should be noted in Eqs. (39) to (42) that the 
inertia factor (1 + y) can always be grouped with the 
time parameter (¢/T,) to give [t/T,(1 + y)]. This 
indicates that if the damping ratio of the system 
f, is kept constant the addition of an inertia component 
y affects only the time of response of the system and 
not the error surges. This is not true of the damping 
component 6, which affects both the response time and 
amplitude of errors of the system. 

Typical solutions of Eq. (39) are given in Fig. 9. 
For an aperiodic control with only an integral com- 
ponent, the error surge following a sudden change in 
engine torque is a full 73 per cent of the deviation with 
Reducing the damping ratio to 0.60 and 


(42) 


no control. 
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Fic. 10. Magnitude of initial error surge. 


making the control oscillatory only reduces the surge 
to 52 per cent. However, the addition of a control 
damping component 6 of magnitude 5.00 cuts the 
initial error surge to 13 per cent while retaining aperio- 
dic response for the control. The advantage of the 
5 component is evident. 

} Variation of the magnitude of initial error surge 
with control parameters is obtained by solving Eq. 
(39) for the time of the error peak and substituting 
this value of time back into Eq. (89) to determine the 
maximum value of n/n,. This is found to be 


3 sin“? V1 - ts? 


(n/N) maz- = [2¢,/(1 + d)Je vi =o (43) 
and this expression is plotted in Fig. 10. Merely 
doubling the available damping by making 6 = 1.00 
reduces the initial error surge as much as reducing the 
damping ratio ¢, to 0.25, and without departing from 
aperiodicity of the control. Eq. (44) further confirms 
the observation made earlier that magnitude of error 
surge after a sudden change in torque is affected only 
by the damping ratio of the system ¢, and the damping 
cofmponent of the control 6 and not by the inertia com- 
ponent of the control y. 


LAG IN THE CONTROL 


The development so far has considered the propeller 
control to be without lag—that is, that an error in 
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propeller speed immediately produced a controlling 
torque proportional to error, its integral, or derivatives. 
However, actual controls have lag that is usually un- 
desirable and which may cause the propeller and control 
to hunt. The incompleteness of the analysis so far is 
demonstrated by Eq. (37), which indicates that the 
system represented by control and propeller is always 
stable and that, when subjected to a transient dis- 
turbance, the system will eventually return to a stabie 
position of zero error, provided only that ¢, > 0. 
Actual controls are known, however, to be capable of 
sustained and self-excited oscillations. The theory 
must then be extended to account for control lag. 

The type of lag in a particular control and its order 
of magnitude will depend on the physical design of the 
control. A few basic types of lag may be recognized. 
It is possible to have control lag proportional both to 
control velocity and acceleration and caused by resili- 
ence or ‘‘spring’’ in the control in the presence of damp- 
ing and inertia, respectively. Constant time lag may 
also be recognized, in which case the response of the 
control lags its error impulse always by a constant time 
interval. 

If the desired operation of the control is represented 
by Eq. (20), introduction of lag proportional to control 
velocity and acceleration can be written 

ae A i 
B=k fndt+...— TA — — B (44) 
w,2 

Notation for the coefficients of Bi and {; is explained 
by noting that an additional degree of freedom has 
been introduced for the control and that, if the control 
is required to wipe out an arbitrarily applied error { in 
blade angle in the case of no speed error m and no 
inertia lag in the control, 8) is reduced according to the 
equation 

B= Be’ (45) 
so that 7’. is the time constant of the control. 

Similarly, if 7, = 0, the control will hunt about zero 


after being constrained to depart 6) from equilibrium 
according to the expression 


B = Bo COS wal (46) 


whence w, is the “natural undamped frequency”’ of the 


control. 
It is further desirable to establish two dimensionless 


lag parameters based on 7°. and w, as 
a = T,/T, (47) 
v = 1/u,T, (48) 
Referring to Eq. (44) the damping ratio of the ffee 
control is noted to-be 
C. = wyl./2 = a/2y. (49) 


Resolving Eqs. (10) and (44) and denoting by p 
the dimensionless operator 





p = T,(d/dt) (50) 
the solution for m written in operational notation is 
n= 


T,(Qo + vQ, (pay? 4 pa + p) 





bY? + PW? + a) +p + y + a) + pl + 6) +(1/4602) 
(51) 





Since Eq. (51) involves the fourth power of )p, in- 
dicating a differential equation of the fourth order, it 
represents motion of a physical system which can easily 
sustain a self-excited oscillation. Applying Routh’s 
rules for stability,’ it is found that Eq. (51) represents 
stable motion if all of the coefficients of the differential 
equation are positive and if 

4&0? > — eon eter a 
(a+ty?)(l+a+y) (1+ 4) — (1+ 6)? 
(: 

Eq. (52) is general for a control of one degree of 
freedom with two orders of lag and both integral and 
derivative components of control. Stability require- 
ments for simpler controls are obtained by letting the 


2) 


~ 


proper parameters go to zero. The case of a simple 
integral control (6 = y = 0) with first- and second- 
order lag is presented in Fig. 11 where stability bound- 
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Fic. 11. Effect of control mass and damping on stability. 


aries are shown. Fig. 11 shows that when only first- 
order lag is present the system is stable for all magni- 
tudes of lag, the only requirement being that {) be 
greater than 0.50. The addition of second-order lag 
for a given amount of first-order lag always makes the 
system less stable, but the converse is not true, Fig. 11 
indicating that an unstable system with inertia lag 
in the control will not always be stabilized by increasing 
the damping in the control. For the latter case, 
stabilization can be obtained only by increasing the 
aerodynamic damping of the propeller ¢) or decreasing 
the inertia lag of the control y. Decreasing the first- 
order lag (control damping) in the presence of second- 
order lag will eventually cause instability in the form 

















of hunting at a frequency approaching the natural 
undamped frequency of the control itself. 

To consider the effect of y and 6 on stability, let 
Eq. (52) be reduced to include only first-order lag. 
Then 


460? > a/(1 + at y)(1 + 4) (53) 


Stability is increased by positive 6 and y. Since 
positive 6 increases the amount of available damping, 
its beneficial effect is to be expected. However, for 
reasons of quick response it is desirable to make y 
negative. The limit of permissible negative y is then 
fixed by Eq. (53). 


MULTISTAGE CONTROL 


The control considered in Eqs. (51) to (53) involved 
only one additional degree of freedom. It is possible 
to have a control consisting of several components in 
series, each with its own degree of freedom. Such 
a control might consist of a centrifugal governor with 
hydraulic power amplification, both governor and oil 
unit having lag. Consider a control consisting of J/ 
units, each subject to first-order lag (7), 72, ... Ts). 
The sequence of control equations is 


Tx -- “= G Ps n dt + oo 
Txt2 + Xo = Cox, 


T3%3 + X3 = CsXe 


T Bi + Bi = CyXy-1 


T,n + n = Qg7T,8 + TyQo (54) 
Noting that 
Cle... Ce = (55) 
the solution for in operational notation is 
T,(Qo + 00.) (tor + 1)... (bawt+ DI (56) 


> = (por + 1)... (baw +1)(p + 1) + (1/ pate?) 


Again making use of Routh’s discriminants, the 
stability requirements for the physical system rep- 
resented by Eq. (54) can be determined for any 
particular value of M. This has been done for several 
cases, and the results are shown in Fig. 12. The case 
of M = 2 is designated as II servo, etc. For M = 2, 
two curves have been plotted, one representing equal 
lag in the two units and the other curve representing 
twice as much lag in one unit as in the other. Note 
that abscissas of Fig. 12 are ay or the sum of the lag 


parameters. That is, 


ri To+... : 
BS eS ee ee ee 


T, 


(57) 
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Then for the case of .\/ = 2, stability of the overall 


system is least for a given total lag a7 when the lag is 
equally distributed Letween the two units of the control. 
A general theorem—stating that when a number of 
servo units each subject to first-order lag are involved 
in a control of the type described, the stability is least 
for a given sum of the lag parameters when the lag is 
equally distributed among the units—is indicated, 
but its proof is beyond the scope of the present paper. 
It suffices to note this fact, and also the additional 
observation from Fig. 12, that for the case of equal 
distribution of the lag among the units for a given 
sum of the lag parameters the stability decreases as the 
number of units over which the lag is distributed in- 
creases. 

The limiting case of an infinite number of servo units, 
each with infinitesimal lag, is obtained by noting that 


b 


limy — ~ [p(ar/M) + 1)” = e**r (58) 


whereupon Eq. (56) becomes 
PT, (Qo + vQ,)e?*? 59) 
= —— cag (oY) 
eT b(p + 1) + (1/40?) 

The physical meaning of the limiting case expressed 
as Eq. (59) is clarified by noting that if the simple 
integral control of Eq. (24) is subject to a constant time 
lag T such that its response at time / is a function of the 
speed error at time ‘—7, then 


(60) 


B, = kin «zy 


Following Minorsky’s’ method of treating constant 
time lag by means of a Taylor’s expansion, Eq. (60) 
becomes 


. ] » 
B= k(n —nT + _ AT? —...)=kme*"r (61) 


Upon combination of Eqs. (61) and (10), Eq. (59) 
results. It is then apparent that the response of WM 
servo units in series, with direct coupling and equal 
first-order lag in each unit, approaches the response 
of a single unit with constant time lag equal to the sum 
of all the lag constants as M approaches infinity. 

Stability of the limit—that is, of the propeller with a 
single unit control having time lag of constant mag- 
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nitude—is determined by assuming that continuous 
oscillation is present in the control. This indicates that 
the equation 


erp(p + 1) + (1/460?) = 0 
has two roots p = *tiw. Substituting into Eq. (62) 
and setting the real and imaginary parts, respectively, 
equal to zero, 


(62) 


(63) 
(64) 


tan apw = 1/w 
wV/1 + w? = 1/4 fo? 


Solving for the minimum value of a, satisfying Eqs. 
(63) and (64) for a given value of {» gives the boundary 
values between stable and unstable motion, and these 
have been plotted in Fig. 12. Evidently, for a given 
lag parameter az, the system is least stable when the 
lag is in the form of a constant time delay in a single 
servo unit. 

Also noted on the curves of Fig. 12 are the simpler 
relations between a7 and {) which the curves approach 
for large values of ap. 


MOTION WITH LAG 


Although stability of the controlled propeller with 
any of the above types of lag present is discovered with 
reasonable facility by means of Routh’s discriminants, it 
is also desired to know something about the character 
of the motion should it prove stable. Further in- 
vestigating the simple integral control with only first- 
order lag, from Eq. (52) the motion is stable if 


4g)? > a/(1 + a) (65) 


The boundary between oscillatory and nonoscillatory 
motion is specified by the requirement for three real 
roots of the differential equation with at least two real 
and equal.’° This is expressed as 


B*(1 — a)? — 2B(1 + a)(2a — 1)(@ — 2) — 27a = 0 
(66) 


where B = 4f;?. 

This boundary is plotted in Fig. 12. 

The speed error after a sudden change in torque is 
given by 


n ei = ; t t 
— = Aye"/To) 4 ¢ nota (A, sin w — + A3 cos] 








Ny T, T> 

(67) 
where 
A; = —(1 — 2re)/[w? + (m1 — 12)?] (68) 
1 ear tha nae) 

Ao = ) 

‘ A w? + (7, — 12)? (m 
2r—1) - 

ha 5 ae (70) 


~ @ + (n — 1)? 


and 7, 72 = iw are roots of the equation in p 
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Pat pil+a)+ pt (1/4602) =0 (71) 


Solutions of Eq. (67) are plotted in Fig. 13. In 
studying the boundary resulting from Eq. (66), two 
questions might arise: first, the penalty for attempting 
to maintain aperiodic motion for an actual control with 
lag; second, the character of the oscillation resulting 
from provision for only sufficient damping to insure 
aperiodic response of an ideal control, although the 
actual control is known to have lag. These two condi- 
tions are shown in Fig. 13. Solid curves show error 
surge of the controlled system with increasing lag as 
aperiodicity is maintained; dotted curves show re- 
sponse as lag is increased for constant ¢>) = 1.00. 


NUMERICAL EXAMPLES 


In order to demonstrate typical values of the parame- 
ters described in the foregoing series of equations, 
the response of the typical engine-propeller combina- 
tion earlier considered has been computed for several 
conditions. Two major conditions are observed: 
flight at 250 m.p.h. at 8,700 ft. and initial take-off 
condition at zero m.p.h. at sea level. Response is cal- 
culated in each case for a control with lag constant 
T, = 1.44 sec. (Cases I and II) and for a control with 
no lag (Cases Ia and IIa). A simple integral control 
is assumed with sensitivity k; = 0.985° per sec. per 
r.p.s. To illustrate the effect of changing control 
sensitivity with blade angle (nonlinear control) the 
take-off condition has also been recalculated for a 
nonlinear control (Case III). Forcing function in 
each case is a sudden 500 ft.lb. change in engine torque. 
For simplicity the torque gradient of the engine has been 


neglected. The following parameters are then ob- 
tained in computation for the cases with control lag. 
Parameter Case I Case II Case III Source 
V 250 m.p.h. Om.p.h. 0 m.p.h. 
Altitude 8,700 ft. sea level sea level ; 
o 0.77 1.00 1.00 std. atm. 
J 1.27 0.00 0.00 J = V/ND 
N 24.1r.ps. 24.1 r.ps. 24.1 r4r.p.s. pee 
n 2.33 2.33 2.33 Eq. (17) 
“ 2.14 0.00 0.00 Eq. (33) 
OC,/OF 1.69 1.69 1.69 en 
8 29° 18° 18° N.A.C.A. 
T.R. 594 

Qn —0.071 —0.022 —0.022 Fig. 4 
qe —0.0023 —0.0018 —0.0018 Fig. 5 
v 0.740 0.260 0.260 Fig. 6 
0p/O0J 15.5 paves shes Fig. 8 
Qn —3.08 —1.24 —1.24 Eq. (18) 
Op —2.40 —2.43 —2.43 Eq. (19) 
cae 0.324 sec. 0.806 sec. 0.806 sec Eq. (13) 
Qo 1.97 1.97 1.97 Eq. (9) 
Ny 0.640 r.p.s. 1.59 rps. 1.59rp.s. Ex. (11) 
a 4.44 1.78 1.78 Eq. (47) 
ky 0.985°/- 0.985°/- 0.066°/- 

sec./r.p.s. sec./r.p.s. sec./r.p.s nase 
fo 1.00 0.40 1.55 Eq. (31) 


The speed error in r.p.m. has been plotted against 
time in seconds in Fig. 14. Comparing Cases I and II, 
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it is noted that the control that produces a mild and 
quickly damped oscillation in flight at 250 m.p.h. still 
allows a violent and continuous hunting at the take-off 
condition. Primary cause of the oscillation is the 
smaller amount of aerodynamic damping available at 
the lower blade angle required for take-off conditions. 
When stability is recovered by making control sensiti- 
vity proportional to blade angle as in the computation 
for Case III, the error permitted by the control becomes 
larger but there is no instability, and for the values 
chosen in Case III control response is aperiodic. 

The harmful effect of control lag is illustrated by 
comparing the error permitted by the controls with lag 
in Cases I and II with the errors permitted under 
exactly the same circumstances by controls of the same 
sensitivity but without lag in Cases Ia and ITa. 


NONLINEAR CONTROLS 


The observation in the preceding computations— 
that with a control of constant sensitivity, response of 
the controlled system may become unstable under 
conditions where blade angle and aerodynamic damp- 
ing are low—leads to the proposal that sensitivity k, 
should vary with blade angle 8. With the aid of Eqs. 
(31), (65), and (66) the boundaries in terms of k; 
and 8 between stable and unstable motion of the 
propeller-engine system for which the curves of Fig. 14 
were computed can be drawn. This has been done in 
Fig. 15. The permissible control sensitivity decreases 
rapidly as blade angle decreases. The stable region is 
larger at sea level than at 8,700 ft. because of the better 
damping in denser air. Increased lag decreases the size 
of the stable region, but for this simple first-order lag 
there is a limit to the unstabilizing effect beyond which 
sufficient aerodynamic damping is present so that no 
amount of first-order lag can produce instability. 

Since the permissible control, sensitivity for low 
blade angles as are required at take-off is almost 
negligible compared with the permissible sensitivity at 
the blade angles of normal flight, the desirability of a 
nonlinear control is evident. Noting the small size 





CONSTANT-SPEED 


PROPELLERS 67 


' 
| —~ CASE 
br Oo mPH 
|? F e SEA LEVE 
“is ff . 2 +18 case 1 
| Ty *!. 44SEC — 
.) 
i / \ wre a vir 
- —CASE 2 \ ‘i Te #144 S80 
 \ 250 MPH \ / K,* 985 
/ 8700 FT / \ / 
\ A-29 \ | 








é 

A 
fe 
X 


| 


a 
| 

| 

l 


gh RS a mat omy 
| 
| 
| 


Computation for actual cases. 





Fic. 14. 


1 














UNSTABLE 
{ REGION 
4 
«© 6¢— 
& / / 
i / fF 
J A 
«4-5 f a 
i 
} y, PA - OSCILLATORY 
© ts Ket (8) 
‘— fr Dg 
La _ somos 
<< —-—- 
P _ (STABLE APERIODIC 
fe} nee PTT TTT a2 P22. WIFI 
iy | | ] | | | 
18 22 26 30 9°34 38 42 
Fic. 15. Stability region for particular propeller-engine com- 


bination. 


of the nonoscillatory region in Fig. 15, it seems also 
desirable to permit some overshooting of the control 
in order to employ sufficient sensitivity to reduce the 
initial error surge as indicated in Fig. 10. 


CONCLUSIONS 


Specific conclusions regarding the particular com- 
ponents of control and their effects are noted in each 
section. In general, it may be said that for propellers 
whose speed is to be controlled by varying the blade 
angle of the propeller the complete system of control 
and propeller will be found most stable at high speed 
and high power, with the possible exception of those 
cases where increase in engine torque with engine 
speed on highly supercharged engines may reduce the 
effective aerodynamic damping and cause hunting. All 
lag in the control is detrimental, leading to increased 
error and hunting, especially at low forward speed of 
the airplane when propeller blade angle islow. Control 
operation is improved by augmenting the aerodynamic 
damping of the propeller by a control component which 
makes change in blade angle proportional to change in 


(Continued on page 70) 











On the Design of the Contraction Cone for a 


Wind Tunnel 
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iy DESICNING the contraction cone for a wind tunnel, 
the usual design condition is that the velocity at 
the end of the cone must be fairly uniform. However, 
if the curvature of the wall along the flow direction 
is too large at certain points, local velocities at these 
points may exceed the uniform velocity at the end 
of the contraction cone. There are then regions of 
adverse pressure gradient and the boundary layer 
may separate from the Furthermore, if the 
velocity at the end of the contraction cone is very high, 
say about 0.9 that of the sound, an additional factor 
namely, the danger of compressibility shock. 


wall. 


enters 
This danger can be avoided by keeping the velocity 
of flow below that of sound in the whole field of flow. 
In the particular case of a contraction cone, the highest 
velocity is reached at the wall of the cone. Therefore, 
if the velocity at the wall is made to increase mono- 
tonicallyt from the beginning of the cone to the end 
of cone, the velocity in the cone will be always less 
than that of sound, provided that the velocity at the 
end of the contraction cone is less than the velocity of 
sound. ‘The pressure along the wall will then be de- 
creasing monotonically. Hence, the danger of boundary 


layer separation is also avoided. 


FORMULATION OF THE PROBLEM 


To actually design such a contraction cone using the 
hydrodynamics of compressible fluids is rather com- 
plicated. However, it is substantiated by both theoretic 
and experimental investigations that the velocity of 
flow of a compressible fluid can be obtained by increas- 
ing the velocity of an incompressible flow with the same 
boundary by a certain factor. This multiplying factor 
is a function of the incompressible velocity itself and 
increases with it.! Hence, if a contraction cone is 
designed for incompressible fluid such that the velocity 
at the wall increases monotonically from the beginning 
to the end, the same cone with compressible flow will 
give a wall velocity that is the product of two mono- 
tonically increasing functions and therefore is itself 
monotone. 

Received September 21, 1942. 

* Research Fellow in Aeronautics. 

7 A monotonic increasing function is one.in which for every 
sequence of increasing values of the argument the corresponding 
values of the function always increase; similarly, a function is 
monotonic decreasing if its values decrease as the argument 


increases. 


To design such a contraction cone for incompressible 
flow is, however, much simpler. One can proceed in 
the following manner: 

At the axis of the cone (Fig. 1), i.e., at ry = 0, assume 

u = fo(x) (ry = 0) (1) 
where wu is the velocity in the x-direction. If v is the 
velocity in the r-direction, by symmetry 

v= O0atr=0 (2 
Values of « and v or the resultant velocity w can be 
found for other values of r by using the following 


fundamental equations of hydrodynamics: 
0v/Ox — Ou/or = 0 (3) 
(0/O0x)(ru) + (0/Or)(rv) = 0 (4) 








Flow in a contraction cone. 


Fic. 1. 


Eq. (3) is the condition for irrotational flow. Eq. 
(4) is the continuity equation. 
By combining Eqs. (2) and (3), one obtains 
ou/Or = 0 atr = 0 (5) 
Hence, u is an even function of ry. Similarly, v is shown 
to be an odd function of 7. Therefore, it is appropriate 


to put 
2 
“om )9r'"s,(x) (6) 
n=0 
o : 
Vv = » gee Gon 41 (X) (4) 
n=0 


By substituting Eqs. (6) and (7) into Eq. (3) and 
equating equal powers of r, 

Lon —1(%) = 2nf2y, (2) (3) 
where the prime indicates derivative with respect to x. 
By substituting Eqs. (6) and (7) into Eq. (4) and 
equating equal powers of r, 




















fn) = = Saee—n(4) (9) 


Eqs. (8) and (9) give the recurrence relation between 
the functions fe, as follows 


a ae" (10) 


Fon(x) = ~ On )? 


Therefore 


— 1)" 
foal) = Seen fo) (11) 





Hence, by substituting these expressions back into 
Eqs. (6) and (7), 





(—1) «abl (2n 
u ‘ 12 
- > For" (x) (12) 
=. (—1)"2nr*"- (2n—1) 
v= 13 
pee i a 
The resultant velocity w can then be calculated as 
w= Vu? + 2? (14) 


The stream lines are lines of constant values of the 
stream function y defined as 


V(x,r) = fC’ru(x,r)dr (15) 


By starting from an assigned monotone velocity dis- 
tribution u = fo(x) at r = 0, the resultant velocity w 
and the stream lines can be calculated by this pro- 
cedure. The shape of the contraction cone is then 
determined by the stream line along which the velocity 
w still varies monotonically, but further out or for a 
further increase in the radius r this condition is no 
longer satisfied. In other words, the shape of the 
contraction cone is determined by the ‘‘last monotone’”’ 
stream line. 

It can be seen that the shape of the contraction cone 
depends upon the function fo(x) chosen at the beginning. 
However, with a plausible assumption of the function, 
it is believed that a satisfactory design for the contrac- 
tion cone can be obtained. 

It is interesting to notice that the solution in the 
form Eq. (12) can be easily identified with the well- 
known Laplace expression for the symmetric potential 
function 


“= L fF folx + ir cos 6)dé (16) 
Tv 


However, for practical numerical calculation the 
series form of Eq. (12) is more convenient. 


SOLUTION OF THE PROBLEM 


To carry out the computation, it is assumed that 
fo(x) = 0.55 + 0.90 - _— e~ */2) dx (17) 
0 


j 
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Assumed velocity u along the axis of the cone’as a func- 
tion of the distance’ on the axis. 
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This distribution of velocity along the axis of the 
tunnel is shown in Fig. 2. The velocity at the end of 
contraction cone is here taken as unity. It is seen 
that the chosen distribution of velocity is quite 





plausible. Furthermore the derivative of fo(x) can 
be easily calculated as follows: 
Put 
&(x) = oe (18) 
rays 
then 
fo" *™ (x) = 0.908 (x) (19) 
But 
a ow ) 
bo” (x epee aan. — (x*/2) 
V/ 25 dx™ ) 
By using the substitutions, x = +/2z, x?/2 = s* 
6” (x) = J eee Sos a &(x)H,,,(z) 
V 24 2!" dx™ gm 


where H,,(z) is the Hermite polynomial. By means of 
the recurrence relation between the Hermite poly- 
nomials, 


H,,(2) = 22H,,-1() — 2(m — 1)H,,-2(2) 


the following recurrence relation between &” (x) is 
immediately obtained 


Be” (x) = — [xb (x) + (m — 1)6™~ (x)] (20) 


The values of §0*(1/+/2r)eo*"dx, (x), (x), 
&'’(x), &’’’(x), 6°” (x) are given in many mathematical 
tables.? The higher derivatives can be calculated by 
using the recurrence relation (20). Then Eq. (19) 
gives all the required functions for calculating u, and 
v is given by Eqs. (12) and (13). For each of the series 
(12) and (13), ten terms are used to give four signi- 
ficant figures. The stream function WV is obtained by 
numerical integration from Eq. (15). 


THE RESULTS 


The stream lines and velocity profiles are plotted 
in Fig. 3. The part of the stream lines where the 
velocity exceeds unity, the asymptotic velocity of the 
cone, is drawn as dotted lines. Therefore, theo- 
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Fic. 3. Stream lines (A B,C, D, E, F) in the meridian x,r 
plane of the contraction cone together with the velocity profiles 
at different sections. 


retically the ‘last monotone” stream line must lie 
inside the line E. However, for the stream line F, 
the velocity at x = 1.2 is 1.032, while at x = 2.0 the 
value is 1.018. This slight excess of velocity will be 
easily removed by wall friction if this stream line is 
used as the wall of the contraction cone. In other 


AERONAUTICAL SCIENCES—FEBRUARY, 


1943 


words, the stream line F can be used as the shape for 
the contraction cone. 

This result has been applied to the design of the 
contraction cone of a large variable density wind 
tunnel. The cross section of the contraction cone is 
not circular but partly circular and partly octagonal. 
Since this modification must tend to increase the local 
velocity at some points of the wall, it is thus thought 
safer to use the stream line D instead of stream line F 
as the basis for designing the contraction cone. Fur- 
thermore, the velocity change from x = 2.0 tox = 
2.8 is negligible. Therefore, the final recommended 
shape of the contraction cone for this particular wind 
tunnel is based upon a boundary as marked in Fig. 3. 
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propeller speed. There seem to be real benefits to be 
derived in the case of integral or “‘rate’’ controls in 
having control sensitivity decrease with blade angle, 
thereby assuring close regulation at high speeds and 
angles and stability at low speeds and angles. 


REFERENCES 


1 Caldwell, F. W., Martin, E., Rhines, T. B., and Baker, C. F., 
The Constant Speed Propeller, Part I, Its Performance, S.A.E. 
Journal, Vol. 40, No. 1, pp. 28-33, January, 1937. 

2 Caldwell, F. W., Martin, E., Rhines, T. B., and Baker, C. F., 
The Constant Speed Propeller, Part II, Development of Its Control, 
S.A.E. Journal, Vol. 40, No. 1, pp. 33-44, January, 1937. 

3 Curtiss Constant Speed Propeller, Aero Digest, Vol. 31, No. 6, 
pp. 60, 62, December, 1937. 


‘Lampton, G. T., Testing of Controllable-Pitch Propeilers, 
A.S.M.E. Trans., Vol. 58, No. 4, pp. 263-266, May, 1936. 

’ Hartman, Edwin P., and Bierman, David, The Aerodynamic 
Characteriscics of Full-Scale Propellers Having 2, 3, and 4 Blades 
of Clark Y and R.A.F. 6 Airfoil Sections, N.A.C.A. T.R. No. 640, 
1938. 

6 Theodorsen, Theodore, Stickle, George W., and Brevoort, 
M. J., Characteristics of Six Propellers, Including the High-Speed 
Range, N.A.C.A. T.R. No. 594, 1937. 

7 Weiss, Herbert K., Constant Speed Control Theory, Journal of 
the Aeronautical Sciences, Vol. 6, No. 4, pp. 147-152, February, 
1939. 

8 Routh, E. J., Advanced Rigid Dynamics, Chapter VI; 
Macmillan Company, London, 1930. 

® Minorsky, Nicholas, Control Problems, Journal of the Frank- 
lin Institute, Vol. 232, No. 5, pp. 451-487, November, 1941; 
No. 6, pp. 519-551, December, 1941. 

10 Dixon, Leonard E., Elementary Theory of Equations, p. 33; 
John Wiley & Sons, New York, 1914. 


The 





Induced Drag of a Twisted Wing 


H. W. SIBERT* 


University of Cincinnati 


INTRODUCTION 


r i ‘HE FORMULA usually given for the induced drag 
coefficient of a wing is (p. 7-28") 


Co, = (1 + o)C,?/2A (1) 


in which C, is the wing lift coefficient, A the aspect 
ratio and o the induced-drag factor. Although valid 
for an untwisted wing (section zero-lift chords all 
parallel), Eq. (1) gives incorrect values for a twisted 
wing. For example, Cp; = 0 in Eq. (1) when C, = 0, 
while the induced-drag coefficient of a twisted wing is 
never zero. 

The correct formula for the induced-drag coefficient 
of a twisted wing can be derived from the theory of 
spanwise air-load distribution given in references 1 and 2. 
In this theory the section lift coefficient C; at any air- 
foil section along the span is expressed as a Fourier 
series, namely (p. 7-29") 


C, = (m,c,/c)(A, sin @ + Ae sin 20 + 
As sin 30 + ...) (2) 


where c is the chord of the airfoil section; c,, the chord 
at the plane of symmetry; m,, the slope of the section 
lift coefficient at the plane of symmetry; A,, a general 
coefficient of the Fourier series; and 0, the arc cosine 
of the distance of the airfoil section from the plane of 
symmetry as a fraction of the semispan. With C, de- 
fined as in Eq. (2), the wing-lift coefficient is (p. 7- 
28") 


Cr = wA UpAy (3) 


where A, is the coefficient of sin 6 in Eq. (2) and u = 
m.c;/4b (p. 7-26"), 6 being the span of the wing. The 
induced-drag factor o is defined by the relation (p. 7- 
28') 

1+o= (A;? + 2A,’ + 3A;? + 4442+ ...)/A;? (4) 
where Aj, As, As, ... are the same A’s appearing in Eq. 


(2). The result of substituting in Eq. (1) the values of 
C, and (1 + o@) given in Eqs. (3) and (4) is 


Co, = tAum?(Ay? + 2Ae? + 343? + 4424+...) (5) 


The correct formula for the induced-drag coefficient 
of a twisted wing can be derived from Eqs. (3), (4), and 


(5). 
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TWISTED WING 


The following analysis treats specifically the case of 
a wing with a defected flap. With a slight change in 
notation, the resulting formulas will also be valid for a 
wing with washout (or washin) at the tips. Let Auy, 
Asy, Asy, ... be the A’s for the untwisted wing (flap 
closed) at an angle of attack of 1 radian from its zero- 
lift position. Let Air, Aor, Asp, ... be the A’s for the 
additional lift due to the flap alone when the ké for the 
flap is one radian. For a twisted wing at an angle of 
attack of x radians from the zero-lift position of the un- 
twisted wing and with the £6 of its deflected flap equal 
to r radians, the corresponding A’s are (xAiy + 1rAip), 
(wAoy + rAop), (xAsy + rAsp), .... The wing-lift 
coefficient and induced-drag coefficient of this twisted 
wing are, respectively (Eqs. (3) and (5)), 


CL = TA Uo(xA iy 4 rA ip) (6) 
wAug?|(xAiy + rAip)? + 2(xAsy + 
rAor)? + 3(xAsy + rAsr)? + ...] (7) 


Let Cry and oy be, respectively, the wing-lift co- 
efficient and induced-drag factor of the untwisted wing 
at an angle of attack of x radians from its zero-lift 
Cry and oy can be found from Eqs. (3) and 
A, by Agu. 


Co; = 


position. 
(4), respectively, by replacing 


Cig = rA UpxA lv (8) 
1 + oy = (Aiy? + 2Aay” + 3Axy? + ...)/Aw* (9) 


From Eq. (9), cy is independent of x and, hence, is 
constant for all values of Cry. 

Let x, be the value of x at which Cp; of the twisted 
wing is a minimum. Cp, is a minimum when 0Cp,/Ox« 
= 0. Hence, from Eqs. (7) and (9), 


each 


Sn = thy (10) 

—(AwAip + 2AoyAor + ...)/(Aw? + 

‘iss... « “ee ee 
SAsyAsp + ...)/Aw*(1 + oy) (11) 


By virtue of Eq. (11), the substitution of this value of 
Xm for x in Eq. (7) gives the following minimum value of 


Coy: 


kn = 


CDim = 9A uo’? K (12) 
; oe = (Air + 2A or’ + ‘+ a + Rm(AiwA tr + 
2AsvAor + ...) = (Air*® + 2A2y* + 


Ra? +...) ~ bth +e) 19) 
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Expanding Eq. (7), 
Co; = wA ug? |x?(A iy? 4 2Aoy? — ea 5 ad 
2Qxr(AwAir + 2AsvAsry + ...) + 


r?(Aip? + 2Aor? + ...)] (14) 


The result of substituting in Eq. (14) the values of 
(Aw? + 2Aey? + ...), (AiwAi + 2AevAor + ..-), 
and (Ajp? + 2Aop? + ...) given in Eqs. (9), (11), and 
(13), respectively, is 


Co, = Au? [(1 + oy) (xAiy — tRmAiw)* +7°Km) (15) 


At this point it will be convenient to introduce the 
quantities Cp;m,’ and Crym’, which are, respectively, the 
minimum Cp, and the corresponding Czy when r = 1, 
that is, when the ké of the flap is 1 radian. From Eqs. 
(12), (10), and (8) 


Com’ = tAug*K,, Clum’ = tAtokmAiy (16) 


The result of substituting in Eq. (15) the values of 
xAyy, RmAiw, and K, given in Eqs. (8) and (16) is 


(17) 


From Eqs. (8), (10), (12), and (16) it is evident that the 
terms r?Cdim’ and rCiym’ in Eq. (17) are, respectively, 
the minimum Cp, and the corresponding Czy when the 
ké of the flap is r radians. This fact checks the validity 
of Eq. (17), since the minimum value of the first term 
of Eq. (17) is zero, which occurs when Cry = rCiym’. 

Let x be the value of x for zero wing lift when the ké 
of the flap is r radians. From Eq. (6) 


x = rko, ko => —Aip/Aw 


Co, = (1 + oy)(Ciy — rCtum’)?/2A + r?Cdim’ 


(18) 


Let Cryo’ be the value of Cry at zero wing lift when r 
= 1—that is, when the ké of the flap is 1 radian. 
From Eqs. (8) and (18) 


Cine’ = wA UpkoA iy (19) 


From Eq. (18) it is evident that rCzy»’ is the value of 
Cry at zero wing lift when the ké of the flap is r radians. 
Hence, if C, is the wing-lift coefficient when the ké 
of the flap is r radians, 


CL = Cry —= rCLypo’ (20) 


The result of substituting in Eq. (17) the value of Cry 
given in Eq. (20) is 
= (1 + oy)[Cr + r(Cryw’ — 
Ciyum'’)|?/xA + r?Cdim’ (21) 


where oy, Cryo’, Ciym’, and Cpi»’ are defined by Eqs. 


(9), (19), and (16). Note that from Eqs. (19) and (16) 
Cryo’ _ Crue’ = TAMA 1y(ko —™ Rm) (22) 


where ky and k,, are defined by Eqs. (18) and (11), 
respectively 
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Example 


Find the formula for Cp, for a wing for which the ké 
of the flap is r radians when the A’s for the untwisted 
wing and for the flap alone are those given on pages 7- 
28 and 7-40, respectively, of reference 1. These A’s 
have been obtained from the procedure outlined in that 


reference. The even A’s were all zero and 


Ay = 0.6777, Asy = —0.0155, Asy = 0.0309 
Az = 0.0009, Asy = 0.0043, 1 + oy = 1.0118 
Aip = 0.4598, Asry = —0.1302, Asr = —0.0495 
Air = 0.0675, Aor = —0.0206 
Ajy? + 2Aoy* + 3Asy? + ... = 0.4647 
Ap? + 2Aor? + 3Asr? + ... = 0.3103 
AyAip + 2AzyvAor + 3AsvAzr + ... = 0.3080 
A = 7.582, u = 0.2801 (pp. 7-24, 7-26) 
ky = —0.4598/0.6777 = —0.67847 (Eq. (18)) 
km = —0.3080/0.4647 = —0.66279 (Eq. (11)) 
K, = 0.3103 — 0.66279 X 0.3080 = 0.1062 
(Eq. (13)) 
Cryo’ — CLum’ = © X 7.582 X 0.2801 X 
0.6777 (—0.01568) = —0.071 (Eq. (22)) 
CdDim’ = m X 7.582 XK 0.2801? XK 0.1062 = 0.1980 
(Eq. (16)) 
(1 + oy)/7A = 1.0118/(4 X 7.582) = 0.0425 
Cp, = 0.0425 (C, — 0.0717)? + 0.1987? (Eq. (21)) 


With this flap deflected 60°, r = 0.1734 (p. 7-32"). 
Hence, for a 60° flap deflection 


Cp, = 0.0425 (C, — 0.0123)? + 0.0060 


WING WITH WASHOUT AT THE TIPs 


Eq. (21) can be applied to a wing with washout (or 
washin) at the tips by deleting r and r? and letting 
CLyo’, CLum’ and Cp,’ refer to the effect of the twist 
alone. Another procedure is to leave Eq. (21) un- 
changed with Ciyo’, Crym’ and Cpjm’ referring to the 
twist alone of a hypothetic wing whose tip chord is 
twisted 1 radian with respect to the root chord. Then, 
yr would be the actual twist in radians of the tip chord 
relative to the root chord. 

The amount of washout (or washin) at the tips in 
most airplane wings is not large and its effect on Eq. 
(21) should be much smaller than the corresponding 
effect of a 60° deflected flap. Hence, for a wing whose 
only twist is a washout (or washin) at the tip, little 
error should result from changing Eq. (21) into 


Co; = (1 + oy) Cr?/rA (23) 
One would hardly be justified, though, in applying 
Eq. (23) to a twisted wing having a flap at its maxi- 
mum deflection. 


REFERENCES 


1 ANC-1(1), Spanwise Air-Load Distribution, Army-Navy- 
Commerce Committee on Aircraft Requirements, April, 1938. 
2'N.A.C.A. T.R. No. 572. 











